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ABSTRACT 


A  method  has  been  proposed  whereby  any  polynomial,  can  be  approxi¬ 
mated  in  an  equiripple  manner  by  a  rational  function.  The  properties  and 
form  of  this  rational  function  are  discussed.  Several  examples  are  used  to 
illustrate  the  theory.  Most  of  these  examples  are  chosen  so  that  the  approx¬ 
imating  rational  function  can  be  identified  as  a  network  function  part;  in 
particular,  the  group  delay  has  been  given  special  emphasis.  The  ideal  group 
delay  vs.  frequency  characteristic  of  filter  is  a  constant.  This  type  of 
group  delay  is  approximated  in  an  equiripple  manner. 

In  addition, a  numerical  scheme  is  proposed  such  that  from  a  given 
crude  equiripple  approximation  a  more  exact  solution  can  be  obtained.  An 
example  used  to  illustrate  this  approach  is  the  problem  of  compensating  the 
non- constant  group  delay  characteristics  of  a  sharp  cutoff  low  pass  filter. 
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1.  HISTCEiCAL  BACKGROUND  AND  STATEMENT  OF  THE  PROBLEM 


1 . 1  Status  of  Equlripple  Rational  Functions  as  Used  in  Filter  Theory 

The  use  of  equlripple  polynomials  (such  polynomials  are  more  commonly 

known  as  Chebyshev  polynomials)  in  filter  theory  has  become  classical  knowl- 
30  3 1 

edge.  ’  Brief ly^  these  polynomials  approximate  zero  in  the  normalized  inter¬ 

val  (-Ij  1)  in  an  equlripple  or  Chebyshev  fashion.  Furthermore,  consistent 
with  the  degree  of  the  polynomial,  the  maximum  number  of  peak  deviations  is 

g 

obtained.  To  attain  this  maximum  number  of  deviations,  all  of  the  polynomial's 
degrees  of  freedom  are  utilized  (the  multiplying  constant  is  not  considered 
a  degree  of  freedom) , 

The  elliptic  function  as  used  in  equlripple  pass  and  stopband 

filter  approximations  is  also  of  a  classical  nature.^  Instead  of  a  polynomial, 

however,  the  elliptic  function  is  rational,  i.e.,  a  ratio  of  two  polynomials. 

Within  the  interval  for  which  zero  is  being  approximated,  this  rational  function 

behaves  like  the  Chebyshev  polynomials.  Outside  of  this  region  a  difference 

occurs  in  that  the  rational  function  also  approximates  infinity  in  an  equiripple 

manner.  Again  the  number  of  maximum  deviations  is  the  largest  number  possible 

21 

for  the  chosen  degree  of  the  rational  function. 

To  obtain  some  freedom  and  more  versatility  in  equiripple  approxi- 

25  5 

mations,  Sharpe  ,  following  closely  the  methods  of  Bernstein.,  developed  a 

rational  function  with  equiripple  properties  within  a  prescribed  interval. 

Furthermore,  the  poles  of  this  function  were  free  to  be  chosen  on  the  imaginary 

frequency  axis  (outside  the  approximating  interval)  and  on  the  real  frequency 

4 

axis  of  the  complex  plane.  Bennett  ,  through  the  use  of  conformal  transformations, 
1.  0 

and  Helman  ,  by  extending  the  definitions  of  the  Chebyshev  polynomials,  suc¬ 
ceeded  in  developing  a  rational  function  with  arbitrary  complex  plane  pole 
locations.  Such  functional  forms  are  of  a  general  type,  the  Chebyshev  poly¬ 
nomials  and  elliptic  rational  functions  being  special  cases.  In  addition  to 

the  above,  many  authors  have  contributed  to  the  area  of  equiripple  functions. 

31 

An  extensive  list  of  some  authors  is  given  in  the  literature. 

1.2  Status  of  Constant  Group  Delay  Approximations 

■A  complete  history  and  a  mention  of  all  who  have  contributed  to  the 
area  of  group  delay  approximations  would  be  a  volume  in  itself.  Thus,  only  a 
few  of  the  most  important  will  be  mentioned. 


2 


The  first  classical  contribution  in  the  area  of  group  delay  approxi- 

28 

mation  was  made  with  the  appearance  of  Thompson's  paper  on  the  maximally 

26 

flat  group  delay.  Later ^  Storch  ^owed  that  the  polynomial  obtained  by 
Thompson  was  a  class  of  Bessel  polynomials,  hence  the  common  name  "Bessel 
filter".  Utilizing  the  great  storehouse  of  knowledge  involving  Bessel  func¬ 
tions,  Storch  gave  this  solution  a  very  elegant  treatment. 

As  the  Bessel  filter  gives  an  excellent  approximation  near  the  origin, 

but  suffers  somewhat  for  higher  frequencies,  a  logical  extension  was  an 

22 

attempt  to  approximate  the  group  delay  in  the  large.  Several  authors,  among 

14  18  10 

the  most  important  of  which  are  Guillemin  ,  Kuh  ,  and  Darlington  ,  have 
succeeded  in  approximating  the  group  delay  in  this  manner  to  give  near  equi- 
ripple  group  delay  characteristics  over  a  larger  frequency  interval.  Briefly, 
their  methods  of  approach  were  (1)  choosing  an  equally  spaced  pole  distri¬ 
bution  parallel  to  the  real  frequency  axis  in  the  complex  plane,  (2)  utilizing 
a  potential  analog  approach  to  determine  the  best  pole  distribution  on  an 
ellipse  in  the  complex  plane,  and  (3)  employing  a  Chebyshev  polynomial  series 
to  approximate  a  linear  phase.  At  this  point  another  approximation  can  be  men¬ 
tioned  which  would  be  classified  as  an  approximation  near  the  origin.  It  is 
the  so-called  Gaussian  magnitude  approximation  due  to  Dishal.^^  The  point  of 
view  used  here  departs  from  that  of  the  previous  authors  in  that  the  phase  is 
attacked  indirectly.  Instead  of  considering  the  phase,  the  magnitude  function 
is  made  to  approximate  the  Gaussian  error  function  by  truncating  its  Taylor 
series  expansion.  The  final  result  is  a  monotonically  decreasing  group  delay 
which  is  relatively  constant  within  the  passband, 

1 6  29 

Following  closely  the  work  of  Helman  ,  Ulbrich  and  Piloty  utilized 
4 

Bennett's  conformal  transformation  to  formulate  a  set  of  non-linear  simultaneous 
equations  that  when  solved  would  give  an  equiripple  group  delay.  The  set  was 
then  solved  using  iterative  techniques  on  a  digital  computer.  The  main  dis¬ 
advantage  to  their  published-results  is  that  the  magnitude  of  the  ripple  is  an 
absolute  number  when,  as  shown  in  tire  next  section,  the  time  domain  error,  or 
dispersion,  is  proportional  (in  a  first  order  approximation)  to  the  fractional 
or  percentage  group  delay  ripple.  Recognizing  this  shortcoming,  the  author 
employed  the  rational  function  developed  in  Chapter  2  to  formulate  a  different 
set  of  non-linear  simultaneous  equations  which  were  also  solved  using  an  itera¬ 
tive  scheme  on  a  digital  computer.  The  author's  work  was  just  completed  on 

2 

equiripple  group  delay,  however,  when  Abele  published  similar  results. 
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1.3  Ideal  Low  Pass  Filters .with  Near  Constant  Group  Delay 

One  objective  of  this  worh  is  to  obtain  several  realizable  transfer 

functions  with  an  equiripple  group  delay  characteristic.  Since  such  systems 

are  used  primarily  for  the  transmission  and  filtering  of  pulses,  a  time  domain 

error  in  the  impulse  response  of  a  filter  due  to  the  group  delay  ripple  would 

be  of  interest.  If  two  reasonably  valid  simplifying  assumptions  are  made, 

an  error  estimate  is  easily  obtained  using  the  method  first  described  by 
32 

Wheeler.  These  assumptions  are  (1)  the  group  delay  ripple  is  cosinu¬ 
soidal,  and  (2)  the  group  delay  is  approximately  constant  over  the  frequency 
spectrum  for  which  the  filter  does  not  give  a  large  attenuation. 

Consider  a  transfer  function  H(jci))  which  has  a  group  delay  t(w)  of 

the  form 

t(co)  =  T  (1  -  (-1)"'  e  cos  -2^  )  (1.1) 

o  CO 

o 

where 

=  average  delay  within  the  filter's  passband 

€  =  the  peak  fractional  group  delay  deviation 

m  =  the  number  of  half  cycles  of  ripple  the  group  delay 

possesses  in  the  interval  (O.w  ) 

'  o 

=  normalizing  constant. 

The  group  delay  lor  such  an  idealization  is  illustrated  in  Figure  1. 

Utilizing  the  definition  of  the  group  delay 


T(w)  m  - 


dB(co) 

dw 


(1.2) 


the  phase  of  H(j(jj)  is  found  to  be 

B(w)  =  -T  Tw  -  (-1)"'  e  -§  sin  1  .  (1.3) 

o  L  nflj  00  J 

o 


23  32 

t’apoulis  has  shown,  using  Wheeler's  method  ,  that  the  impulse  response  i(t) 
of  the  transfer  function  H(joo)  with  the  phase  characteristic  of  (1.3)  can  be 
put  in  the  form: 


J.  (6)  h  (t  -  —  ) 
k  00 


i(t) 


(1.4) 


4 


where 


h(t)  is  the  inverse  of  |  H( jcj)  |€ 


-JWT 


J,  is  the  Bessel  function  of  the  first  kind,  and 
k  ’ 

er  w 

c  o  o  , 

6  -^r-  (-1)  • 


If  6  is  very  small  compared  to  one,  (1.4)  reduces  to 


i(t)  ~  h(t)  +  -I  h(t  -  —  )  -  I  h(t  +  —  ) 

2  (jl)  2  CO 

o  o 


(1.5) 


For  a  lumped-element  system  the  phase  of  the  transfer  function 
approaches  n7I'/2  as  oj  tends  to  infinity,  where  n  is  the  number  of  poles  less 
the  number  of  finite  zeros.  Hence,  the  area  under  the  group  delay  curve 
from  u  =  0  to  oj  =  “o  is  nir/2.  If  the  fractional  part  rj  of  this  area  follows 
the  boundary  of  the  area  from  w  =  0  to  cj  =  of  Figure  1,  then 


T 

O 
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and  (1.5)  can  be  written  as 


i(t)  =  h(t)  +  (-I)"*  ^  h  (t 


V 


A  typical  impulse  response  will  then  appear  as  shown  in  Figure  2. 

From  the  foregoing  analysis,  several  conclusions  can  be  drawn.  Some 
of  the  more  important  are: 

1.  The  output  pulse  is  dispersed  in  the  sense  that  the  time  bandwidth 
product  is  increased  because  of  the  non-constant  group  delay. 

2.  The  leading  "echo"  appears  to  give  an  anticipatory  feature  to 
the  filter.  In  a  physical  system,  however,  this  echo  would 

in  effect  cancel  with  the  main  pulse  to  prevent  such  a  phenomenon. 
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Figure  2.  Impulse  response  of  a  filter  with  the  idealized  group  delay  of 
Figure  1 
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3.  The  magnitude  of  the  two  main  distortion  terms,  centered  at 

±  2j^iri/n  is  proportional  to  the  fractional  ripple  e. 

4.  The  distortion  or  error  terms  are  proportional  to  n/m,  i.e., 
the  ratio  of  the  number  of  T^/2  multiples  of  phase  at  to  =  ®° 
to  the  number  of  half-cycles  of  ripple. 

5.  The  dispersion  terms  are  proportional  to  the  constant  jj. 

1.4  Statement  of  the  Fh’oblem 

The  Chebyshev  rational  function  of  Helman  and  Bennett  is  a  very 
practical  and  useful  form  as  applied  to  network  problems.  However,  general 
and  useful  as  it  is,  it  does  possess  some  basic  limitations.  Two  of  these 
are: 

1.  If  a  polynomial  more  general  than  a  constant  is  to  be  approxi¬ 
mated  in  the  Chebyshev  sense  a  serious  difficulty  is  encountered. 
For  example:  using  Bennett's  conformal  transformation,  suppose 

a  linear  term  is  to  be  approximated  with  a  rational  function 
of  a  chosen  degree.  The  end  result  demands  that  the  poles  have 
to  be  chosen  in  such  a  manner  that  a  zero  will  appear  in  the 
proper  location  to  nullify  the  effect  of  a  particular  pole  at 
the  same  location.  Because  of  its  form,  this  would  indeed  be 
a  difficult  task. 

2.  Some  problems  require  a  rational  function  approximation  of  a 
polynomial  such  that  the  number  of  ripples  is  less  than  the 
maximum  possible  number.  Again  the  currently  available  Chebyshev 
rational  function  is  not  adaptable  to  such  a  situation. 

The  problem  is  thus  to  develop  a  closed  form  rational  function  which 
is  general  enough  to  solve  many  of  the  problems  solvable  with  the  existing 
Chebyshev  rational  function  and,  in  addition,  be  sufficiently  flexible  to 
overcome  the  above  stated  difficulties.  Once  such  a  function  is  developed 
there  exist  many  practical  network  problems  to  which  it  can  be  applied;  some 
solved,  and  some  unsolved.  Thus,,  after  the  main  problem  has  been  studied, 
several  of  the  solved  and  unsolved  varieties  can  be  investigated  with  the 
purpose  of  exploiting  the  potentials  of  another  tool  in  today's  existing  body 
of  knowledge.  Special  emphasis  will  be  given  to  group  delay  problems,  as  it 
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has  been  only  recently  that  equiripple  solutions  to  this  problem  have  appeared  ' 
and  another  method  of  attack  would  be  of  value. 

As  general  as  the  rational  function  is  which  overcomes  the  two  stated 
difficulties  of  the  Bennett  and  Helman  forms,  it  is  not  difficult  to  think  of 
equiripple  or  weighted  equiripple  problems  which  are  not  easily  adapted  to  fit 
any  equiripple  rational  function.  With  this  in  mind,  a  numerical  scheme  is 
suggested  and  illustrated  to  solve  some  of  the  existing  group  delay  problems 
which  may  be  encountered. 
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2.  THE  DEVELOPMENT  OF  A  PROCEDURE  FOR  APPROXIMATING  POLYNOMIALS 
WITH  rational  functions  TO  GIVE  AN  EQUIRIPPLE  ERROR 

Many  problems  arise  wherein  a  rational  function  is  desired  which 
approximates  a  polynomial.  Often  this  polynomial  is  of  a  low  degree;  indeed, 
most  practical  problems  demand  that  it  be  nothing  more  than  a  constant.  The 
problem  undertaken  in  this  chapter  is  thus  the  determination  of  some  proper¬ 
ties  and  how  to  construct  si.ch  a  rational  function. 

2.1  The  General  Case 

The  rational  function 


R 

nm 


(X) 


n 


z 


b .  X 
1 


i=o 


(2.1) 


of  given  numerator  and  denominator  degrees  n  and  m  respectively  is  to  be 
found  which  approximates  the  given  polynomial. 


i 

P^(x)  =  ^  (2.2) 

i=o 

of  degree  Z  in  an  equiripple  manner  in  the  interval  a  <  x  <  S.  The  constants 

a.,  b. ,  and  c.  are  all  real,  and,  without  loss  of  generality,  b  is  assumed  to 
i’  1  ’  ’  ’  m 

be  unity.  If  R  is  to  approximate  P,,  in  [a, B]  with  an  equiripple  error  of  peak 
nm  jt 

deviation  £  then 


P-e<R  <P^+t  lor  a  <  X  <  3. 
Z  -  nm  —  Z  -  - 


10 


The  number  of  times  the  equality  is  satisfied  as  x  takes  on  all  values  between 
CL  and  B  will  be  called  n^ .  An  upper  bound  on  n^  will  be  determined  subsequently. 

As  an  intermediate  step  in  determining  the  numerator  A  (x)  and  denom- 

n 

inator  B  (x)  or  the  rational  function  R  (x) ,  the  error  functions  E  (x)  and 
m  nm  ■’  + 

E_ (x)  will  be  defined.  From  these,  some  useful  properties  will  be  deduced. 


The  one  exception  occurs  when  n  =  m  +  X  (n  ^  0)  and  a  =  c. .  For  this  com- 

n  At 

bination,  the  degree  of  and  E_  is  k  -  1.  If  n  =  m  +  ^  with  n  =  0  the  prob¬ 
lem  is  trivial  since  only  constants  are  involved.  Thus,  it  can  be  concluded 
that  the  number  of  zeros  of  E^  and  E  is  k  (for  the  above  noted  exception  k  -  1). 
For  R  to  approximate  P,  in  an  equiripple  manner  with  the  largest  possible 
number  of  maximum  deviations  of  i  €  all  zeros  of  E^  and  E  are  required  to  be 
real.  Furthermore,  because  of  the  tangency  at  the  internal  zeros  these  will 
be  double.  The  two  zeros  at  the  ends  of  the  approximation  interval  wfil  be 
simple . 
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The  maximum  number,  n^,  of  i  €  deviations  is  now  easily  determined. 
Knowing  that  the  total  number  of  zeros  involved  in  and  E_  is  2k  and  all  but 

two  of  these  are  double,  the  number  of  extrema  points  is 


n 

e 


+  2 


k  +  1  , 


Unless  stated  otherwise  in  the  following,  the  k  zeros  of  E^  and  E 
will  be  assumed  to  be  (1)  real,  (2)  all  double  except  the  two  at  each  end  of 
the  interval,  and  (3)  interlaced,  i.e.,  in  going  from  x=atox=Ba  zero 
of  E^(E_)  will  be  followed  by  a  zero  of  E_(E^),  etc. 

Making  the  definitions 


numerator 

of 

E 

+ 

=  A  - 
n 

■ 

+  e)B 

] 

numerator 

of 

E 

=  Q_ 

=  A  - 
n 

-  6)B, 

it  is  possible  to  distinguish  five  different  cases  which  will  exist  depending 
upon  the  degrees  i,  m,  and  n.  These  cases  are  listed  as  follows; 


Case  la;  n>i+m  ;  i  >  1  ;  k=n 


^  n  ■  ^i^ 


i=l 

k 


^  n  -  "'i^ 


i=l 


Case  IB:  n=f+m  ;  i>l  ;  k=n 


k  . 


^  n  ■  «i^ 

i=l 

k 

-  <•„  -  'i*  n 


i=l 


(2.6) 


(2.7) 


(x  -  h^) 
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Case  II;  n>m  ;  X=0  ;  k=n 


k 

\  n 

i=I 

k 

u 

i=l 

Case  III;  n=in  ;  JL  =  0  ;  k  =  n 

k 

n  -  e.) 

i=l 

k 

Q  =  (a  -  c  +  £)  TT  (x  -  h  ) 
-  no  J-J-  1 

1=1 

Case  IV;  n<m  ;  Z  =  0  ;  k  =  m 

k 

=  -  (c^  +  €)  n  -  Si) 

1=1 

k 

Q  =  -  (c  -  €)  n  (X  -  h ) 

—  O  1 

1=1 

Case  V;  n<m  +  ;^  ;  Z  >  1  ;  k=m  +  £ 

k 

%  =  -  ^zU  ■  ^1^ 

1=1 

k 

=  -■  n 

1=1 


(2.8) 


(2,9) 


(2,10) 


(2,11) 
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Two  important  differences  to  note  in  the  above  cases  are  (1)  whether  or  not 
the  constant  multiplying  factors  of  Q  and  are  the  same  and  (2)  the  value 
of  k. 


Knowing  Q  and  in  terms  of  the  multiplying  constant  and  the  points 

of  maximum  deviations,  it  is  a  simple  matter  to  solve  for  A  and  B  from  (2.5), 

n  m 

The  results  are 


A 

n 


(p^  +  €)  Q_  -  (P^  -  e)  Q 
__ 


(2.12) 


B 

m 


2€ 


Hence,  R  is 
’  nm 

R  ,  ^  ^  Q_  -  (P,  -  Q, 

nm  B  Q  -  Q 

m  -  + 


(2.13) 


which  can  also  be  rearranged  to  fit  into  the  forms 


Q  +  Q, 


R  =  Pfl  +  €  — 
nm  X  Q 


=  P«  +  e 


R 

nm 


=  P- 


€  +  26 


(2.14a) 


(2.14b) 


(2.14c) 


The  number  of  arbitrary  constants  for  a  given  case  can  now  be  deter¬ 
mined.  The  solution  will  be  given  in  detail  for  Case  la.  The  results  are 
tabulated  in  Table  I  for  the  remaining  cases. 

The  defining  relationships  of  Case  la  are  n>X  +m  ;  X>1,  From 
these  it  is  evident  that  k  =  n,  and  n^  =  n  +  1.  The  following  information  is 
obtained  by  an  inspection  of  (2,1),  (2.2),  (2.6),  and  (2.13). 
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n  £  degree  ofQ  -  Q  =n-l 

n  H  degree  of  Q  (P.  -£)  =n+X-  1 

4  4" 

n_  =  total  number  of  constants 

5 

=1  (e)  +  n  (number  of  error  extrema) 
e 

+  f  +  1  (number  of  constants  in  P^) 

4-  n  (n  =1  for  case  2,  zero  otherwise) 

=  n  +  fr  +3 

n„  =  number  of  numerator  coefficients  to  be  zero  so  that  n  is  the 

6 

desired  integer 
=  n^  -  n 
i  -  1 

n^  =  number  of  denominator  coefficients  to  be  zero  so  that  m  is  the 
desired  integer 
=  ng  -  m 
=  n  -  1  -  m 

n„  =  number  of  constraints  which  have  to  be  satisfied 

O 

=  +  n„ 

6  7 

=  i  +  n  -  m  -2 

Ug  =  number  of  arbitrary  constants 

=  ^9  -  ”8 
=  m  +  5 

A  word  of  caution  should  be  inserted  at  this  point.  The  arbitrary 
constants  are  not  unlimited  in  value.  Rather^  there  are  definite  bounds  placed 
upon  their  range,  this  range  being  defined  by  the  particular  problem  being 
solved. 

An  example  will  now  be  given  to  illustrate  the  foregoing  theory.  A 
rational  equiripple  approximation  to  the  polynomial  P^  =  x  is  desired  in  the 
interval  (1,  10).  The  error  extrema,  €,  is  to  be  1/2,  and  the  degree  of  the 
numerator  and  denominator  is  to  be  4.  A  solution  is  as  follows: 
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Since  X  =  1  and  m  =  n  =  4,  the  inequality  n  <  £  +  m  is  satisfied 
and  the  problem  fits  Case  V,  hence, 

(a)  the  maximum  number  of  extreme  points  =  6 

(b)  the  number  of  arbitrary  constants  =  8 

(c)  the  number  of  constraints  =  1. 

Five  of  the  arbitrary  constants  have  been  specified  (£,  c^,  c^^,  a,  B) .  Since 
the  number  of  extrema  —  of  which  two  are  specified  and  three  are  arbitrary  — 
is  6,  one  extreme  point  will  be  defined  when  the  others  are  chosen.  The 
approximation  will  be  similar  to  the  illustration  of  Figure  4, 

From  (2.13),  the  general  form  of  R  is 
’  nm 


Q_  (P^  +  6)  -  (P^  -  £) 

*^44  ""  Q_  - 

Using  the  notation  shown  in  Figure  4 

=  -  (x  -  1)  (x  -  x^)^  (x  -  x^)^ 

2  2 

Q_  =  -  (x  -  10)  (x  -  Xg)  (x  -  x^) 


(2.15) 


(2.16) 


The  requirement  that  n  =  4  places  a  restriction  on  the  choice  of  the  set 

5 

(x^,  x^) .  Equating  the  coefficient  of  x  to  zero,  or  setting 


^4  - 


^2  - 


x^  -  4  =  0 


(2.17) 


satisfies  the  demand. 

An  arbitrary  choice  of  x^,  x^,  and  x^  is  3,  4,  and  8  respectively. 
This  choice  makes  x^  =  11,  but  this  is  outside  the  approximation  range  and  so 
cannot  be  used.  Similarly,  a  choice  of  1,  7,  and  8  for  x^,  x^,  and  x^  respec¬ 
tively  requires  x^  =  6,  Again  this  cannot  be  used  because  must  be  greater 
than  x^  and  less  than  10.  A  choice  of  3,  5,  and  7,  however,  gives  a  realiz¬ 
able  x^,  namely  9,  With  the  last  set  of  x^'s  the  rational  function  reduces  to 


R 


44 


1 .5x 


4 


4 

X  - 


34x^  +  337x^ 
28x^  +  294^  - 


1578x  +  3217.5 


(2.18) 


1356X  +  2385 
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Before  continuing  on  to  another  topic,  one  more  detail  needs  to  be 
investigated.  This  detail  concerns  the  possibility  of  the  existance  of  another 
rational  function  which  (1)  passes  through  ^  £  at  the  ends  of  the 
approximating  interval  (a,  6),  (2)  is  still  tangent  to  and  passes  through 
P.  i  £  at  the  internal  points  x  ,  x  ,  ...  x  ,  and  ,(3):  yet,  doesn't 

a  2 

intersect  P^  within  (a,  6) .  Seemingly,  such  a  rational  function  (Figure  5) 

is  easily  sketched.  From  such  a  sketch,  however,  it  is  observed  that  outside 

the  interval  (a,  6)  the  function  R  must  satisfy 
'  nm 


1^ 


<  £ 


(2.19) 


or  R  will  intersect  with  P .  i  £ .  Such  an  intersection  is  impossible, 

however,  because  the  degree  of  the  numerator  of  is  not  large  enough  to 

permit  any  more  intersections.  As  £,  a,  and  6  are  arbitrary,  it  is  concluded 

that  the  inequality  (2.19),  cannot  hold..  Hence,  R  does  not  exist. 

'  nm 

2.2  Even  Rationals  Approximating  Even  Polynomials 

The  case  under  consideration  is  a  special,  but  important,  case  of 
Section  1.  Because  both  R  and  P.  are  even,  only  one  half  of  (-S,  6)  need 
be  considered.  The  approximation  is  carried  out  as  though  P^  was  being  approxi¬ 
mated  in  (0,  3).  Then  at  every  zero  of  E  and  E  the  factor  (x  -  x.)  is  replaced 
2  2  +  -  1 
by  (x  -  x^  ) ,  Several  examples  of  this  type  of  approximation  are  given  in 

Chapter  3. 


2,3  Odd  Rationals  Approximating  Odd  Polynomials 


Again  this  is  a  special  case  of  Section  1.  Under  the  assumption 

that  R  and  P.  are  odd,  the  following  conclusions  can  be  drawn: 
nm  X  ’ 


that 


a.  From  the  arbitrary  odd  approximation  shown  in  Figure  6  it  is  noted 
K_  Q_  (x)  =  -K^  (-X) 


where  and  K  are  multiplying  constants. 

b.  Since  P^  is  odd,  X  /  0;  hence  from  Table  1  it  is  seen  that  the  • 
multiplying  constants  of  Q_  and  are  equal.  Thus 

K  =  K  . 

+ 
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X 


Figure  5.  The  rational  function  R  and  R  which  both  go  through 

at  X  =  o,  8  and  are  tangential  ?o  ^  ,  and  pass  through 

at  the  Internal  points  x, ,  x„,  x 

1'  2'  ne-2 


H*  H“ 
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c.  Finally,  it  is  observed  that 

p^(x)  +  e  =  -  [p^(-x)  -  e]  . 

Putting  all  three  conclusions  into  (2.13) 

Q^(-x)  [P,(-x)  -  €]  -  Q^(x)  [P  (x)  -  e] 

_  _± _ ± _ j: _ ^ _ 

nm  -Q^C-x)  -  Q^(x) 

Odd  {Q^(x)  [P_g(x)  -  £]} 

"  [Q^(x)] 


(2.20) 


where  Odd  and  Ev  are  operators  which  are  defined  to  take  the  odd  and  even 
parts,  respectively,  of  what  follows. 

Examples  of  this  type  approximation  are  given  in  Chapter  3. 

2.4  Observations  and  Summary 

In  this  chapter  a  rational  function  has  been  developed  which  will 
approximate  any  polynomial  in  an  equiripple  manner.  The  method  has  several 
desirable  features.  Some  of  these  are: 

a.  There  is  a  large  amount  of  freedom  available  in  the  form  of 
arbitrary  constants,  for  fitting  the  rational  function  to  a 
wide  range  or  problems. 

b.  The  rational  function  need  not  he  restricted  to  approximating 
only  constants  as  does  the  rational  function  of  Helman  and 
Bennett . 

c.  Though  not  mentioned  explicitly,  this  rational  function  need 
not  approximate  P^  with  the  maximum  number  of  ripples.  If 
fewer  ripples  are  desired,  instead  of  using  the  factor  (x  -  x^) 
in  E^  where  x^  is  positive  and  real,  the  term  x^  can  be  made 
negative.  Indeed,  quadratic  factors  with  complex  roots  can 

be  used,  if  desired,  to  decrease  the  number  of  ripples.  As 
the  constants  within  these  factors  are  arbitrary,  the  rational 
function  is  given  a  more  arbitrary  character.  An  example  of 


this  type  is  given  in  Chapter  3 


I 


Sbme  limitations  on  R  as  developed  in  this  chapter  are; 

nm 

a.  The  denominator  is  the  difference  of  two  polynomials  and  is 
not  in  factored  form. 

b.  In  developing  a  rational  function  R  which  approximates  a  high 

nm 

order  polynomial  or  has  several  zeros  at  infinity,  several 
noh-linear  simultaneous  equations  will  have  to  be  solved. 

The  above  limitations,  however,  are  overcome  to  a  large  extent  by 
the  ever  increasing  availability  of  high  speed  digital  computers. 


3.  APPLICATIONS  OF  THE  EQUIRIPPLE  ERROR  RATIONAL 
FUNCTION  OF  CHAPTER  2  TO  NETWCEK  PROBLEMS  I 
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The  well  known  theorem  of  Weierstrass  says  that  any  function  con¬ 
tinuous  in  the  interval  (a,  b)  may  be  approximated  uniformly  by  polynomials 
9 

in  this  interval.  In  general  this  requires  a  polynomial  of  infinite  degree. 
Fortunately,  there  exists  a  large  class  of  problems  which  do  not  require  an 
approximation  of  such  high  degree.  Indeed,  in  many  problems  the  desired 
function  is  itself  a  low  order  polynomial.  Such  is  the  case  with  many  of  the 
network  function  parts. 

A  network  function  is  a  rational  function  of  the  complex  variable 
s,  which  may  in  turn  be  broken  up  into  several  parts,  many  of  which  are  them¬ 
selves  rational  in  the  real  variable  ou.  Among  these  are  the  real  part,  the 
imaginary  part,  the  magnitude  squared  function,  and  the  group  delay  function. 
This  chapter  will  thus  be  devoted  to  finding  a  rational  which  approximates 
a  low  order  polynomial  and  yet  can  be  identified  with  one  of  the  above  network 

function  parts.  This  in  turn  may  be  physically  realized  by  conventional 

30 

techniques  if  desired. 


3 . 1  Real  Part  Approximation 

The  real  part  of  a  network  function  H(s)  is  a  rational  function  of 
the  real  variable  od.  That  this  is  the  case  can  be  shown  very  easily  by  con¬ 
sidering  the  network  function. 


H(s) 


+  Nj^ 


(3.1) 


where  and  N^  are  the  even  and  odd  parts,  respectively,  of  a  Hurwitz 
polynomial  in  the  complex  variable  s  =  O  +  juu.  Similarly,  and  N^^  are  the 
even  and  odd  parts  of  any  polynomial.  From  (3^1)  the  even  part 


“l“2  -  ''l''2 
2  2 
“2  -  ^2 


Ev  H(s) 


(3.2) 
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is  easily  found.  By  inspection  this  is  also  a  rational  function.  On  the 
real  frequency  axis  s  =  ju),  and  the  even  part  becomes  the  real  part 


Re  H(ju)) 


“i“2  -  ""iNa 
2  2 
“2  -  ^2 


S=jU) 


(3.3) 


which  is  rational  in  the  real  variable  u). 

If  a  rational  function  is  to  approximate  a  polynomial  and  yet  be 

identified  as  a  "positive  real"  real  part  it  must  not  only  approximate  the 

30 

polynomial,  but  satisfy  the  following  properties; 

a.  possess  no  poles  on  the  imaginary  axis 

b.  the  poles  and  zeros  must  have  quadrantal  symmetry 

c. ,if  the  denominator  of  H(s)iS  of  degree  t,  then  the  denominator 

of  Ev  H(s)  is  of  degree  2t.  The  degree,  2r,  of  the  numerator  of 
Ev  H(s)  must  satisfy  the  inequality 

0  <  2r  <  2t 

d.  imaginary  axis  zeros  of  Ev  H(s)  must  be  of  even  multiplicity. 

If  H(s)  is  not  positive  real  d  and  a  can  be  relaxed. 

Several  examples  will  now  be  given  wherein  a  low  order  polynomial 
is  approximated  by  the  rational  function  of  Chapter  2,  and  yet  is  manipulated 
to  satisfy  some  or  all  of  the  above  positive  real  properties. 

Example  (3.1) 

A  rational  function  is  to  be  obtained  which  satisfies  all  of  the 
above  positive  real  properties  and  in  addition  is  to 

a.  approximate  the  degenerate  polynomial  "unity"  with  a  tolerance 
of  unity 

b.  have  a  maximum  number  of  extreme  points  which  are  separated  by 
equal  frequency  intervals 

c.  have  a  denominator  of  degree  2N 

d.  approach  the  constant  value  of  4  for  high  frequencies. 

An  illustration  of  the  desired  function  is  shown  in  Figure  7  for  N  even  and 
N  odd. 
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Solution: 


It  is  first  observed  that  X  =  0,  and  m  =  n  so  that  Case  III  of 

Chapter  2  applies  to  this  example.  Furthermore,  R  and  P.  are  both  even 

nm  Jl, 

so  the  corresponding  simplifications  can  be  used.  For  the  half-interval 

(0,  N)  Nii 

2 

=  (a^j^  -  1  -  €)  (uj  -  N)(u  Y[  [cD  -  (N  -  2i)]^ 

i-1 


N/2 

"  ^^2N  ■^■*■^^11  [‘U  -  <N  -  2i  +  1)]^ 

i=l 


>  ;  N  even 


(3.4) 


N-1 

2 


%  "  ^^2N  -  1  -  N)  n  -  (N  -  2i)]‘ 


i=l 


M-X 

2 


S  ;  N  odd. 


=  (a^j^  -  1  +  €)uj  [uJ  -  (N  -  2i  +  1)]‘ 


i=l 


The  constant  a  can  be  used  to  satisfy  requirement  (d) ;  then  using  (2.13)  with 
6=1  the  final  result 


e  N 


N/2  2 

4  n  [u)^  -  (N  -  2i  +  1)^] 
i=l 

N-2 

N/2  2  2 

n  [0)"  -  (N  -  2i  .  1)2]  -  -  n2)  n  [0^2 

i=l  i=l 


(3.5) 

-  ;  N  even 

2  2 

(N  -  2i)2] 


2 
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4  [u)^  -  (N  -  21  +  1)^] 


R  H^(uu^)  =  - ^ 

e  N  N-1 


;  N  odd 


[tjo  -  (N  -  2i  +  1)  ]  -  (o)^  -  N  )  -  (N  -  2i)  ] 


valid  in  [-N^N]  is  obtained.  With  N  =  1,  2^  3^  and  4^  the  above  reduces  to: 


R  H  (uu^)  =  - 

el  2 

U)  +  1 


■>,“2  <“'> '  'K ; 

CO 


R  H  (cuS  -  40)^  (0^^  -  4)^ 

»  "  -  5»n  13a,"  t  9 


“,«4  <“'>  ■  8"  “‘"e  '  '  "-s - 

tu  -  16(u  +  92(1)  -  104(1)  +  162 


By  inspection,  (3.6)  satisfies  the  properties  of  the  real  part  of  a  positive 
real  driving  point  function. 

Example  (3.2) 


A  rational  function,  G  ,  is  to  be  obtained  with  the  following  con¬ 


straints: 


a.  approximates  the  degenerate  polynomial  zero  in  an  equiripple 


b.  lim  R  ,G„  =  0 
(1)-Kxi  e  N 

c.  R  G„(n^)  =  (-1)^"  for  n  =  0,  ±  1,  ±  2,  .  .  ±  (N-1) 

e  N  f  f  *  f 


d.  R  G  =  0 

e  n  V22y 

e.  has  2N-1  extreme  points. 
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The  rational  function  to  be  determined  is  sketched  for  positive  frequencies 
in  Figure  8. 

Solution: 

This  problem  appears  to  fit  Case  IV,  however,  two  differences  are 
noted;  (1)  the  approximation  interval  extends  from  minus  to  plus  infinity,  and 
(2)  the  number  of  extreme  points  is  2N-1  which  is  less  than  the  maximum  2N+1. 

To  modify  Case  IV  such  that  is  meets  the  stated  requirements, 
imaginary  zeros  are  placed  at  ui  =  1  ja  in  (Q_)  for  N  even  (odd),  thus 
giving: 

N-2 

^  2 

^  2.2  2^  r-r  .  2  2. 

=  -0)  (a'  +  ^  )  J_J_  -  4i  ) 

i=l 

N/2  2 

Q_  =  n  -  (2i  - 1)^] 

i=l 


>  ;  N  even 


N-1 

2 

„  2  rr  ,  2  _2  ‘ 

=  -UJ  (u)  -  4i  ) 

i=l 


N-1 

2 

Q_  =  (u)^  +  a^)  [u)^  -  (2i  -  1)^] 

i=l 


;  N  odd. 


(3.7) 


Substituting  these  values  in  (2,13) 


N/2 


^  [u)  -  (2i  -  1) 


N-2 

2  ^2 
2]  _  .  2  2.  -rr  .  2  ..2. 


e  N 


i=l 


(uj  +  a  )  (u)  -  4i  ) 

i=l 


N-2 

2 


N  even 


(3.8) 


N/2  2 

[u)^  -  (2i  -  1)^]  +  U)^  (u)^  +  a^)  PJ  (OD^  -  4i^) 

i=l  i=l 


2 
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N-1 

2 


N-1 

2 


(»"  +  n  .  (21  .  1,2]  .  („2  .  4i2,' 


R  G  = 
e  N 


i=l 


1=1 


N-1 

2 


N-1 

2 


;  N  odd.  (3;8) 


+  a^)  n  -  (21  -  1)^]  +  U)^  n  (“^^  -  4i^) 


1=1 


1=1 


As  a  Is  arbitrary  It  can  be  specified  so  that  the  numerator  of  (3.8)  Is  zero 

at  u)  =  1/2.  After  several  algebraic  steps,  the  value  of  N(N-l)  Is  obtained 

^  2 
lor  a  . 

Later  In  this  chapter  will  be  used  for  another  example,  there¬ 

fore,  it  is  tabulated  for  N  =  2  through  6  as  follows: 


R  G 
e 


2  ~ 


-4uj  +  1 
4 

2oo  +  1 


\°3  = 


4  2 

12tu  -  27co  +6 

6  4  2 

2uu  -  4u)  +  5u)  +6 


^  ^  -24uj^  +  198u)'^  -  372(ju^  +  81 

®  ^  2u)^  -  160)^  +  38(u”^  +  120)^  +  81 

40uj^  -  810u)^  +  47400^*^  -  7615u)^  +  1620 
2u).^°  -  40(1)®  +  246u)®  -  380u)^  +  577u)^  +  1620 


(3.9) 


-60(1)^^  +  2415u)®  -  31860(11®  +  155535u)”^  -  239430u)^  +  50625 

1.2  10  8  0  4  2 

2(1)  -  80(1)  +  1071(1)  -  5300(1)  +  10127(i)  +  6330u;  +  50625 


It  is  apparent  that  the  above  functions  satisfy  the  properties  of 
a  real  part  function;  however,  it  is  not  a  positive  real  one. 

Example  (3.3) 

A  rational  function  R  K„  is  to  be  found  which  has  the  properties  of 

e  N 

a  real  part  function  and  yet  approximates  cos  n  (i)  within  the  interval  (-1,  1) . 
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Solution: 

The  function  developed  in  Example  (3,1)  approximates  cos  rr  U)  in 

the  required  interval;  however,  a  better  approximation  can  be  obtained  if 

the  rational  function  is  zero  not  only  at  u)  =  -  1/2,  but  also  at  u)  =  -  3/2. 

To  satisfy  this  last  added  condition  requires  another  arbitrary  constant. 

A  rational  function  is  thus  chosen  with  the  maximum  number  of  extreme  points 

all  except  the  last  of  which  are  separated  by  a  unit  increment.  The  last 

extreme  point  and  the  value  of  at  infinity  will  then  be  used  to  place 

two  zero  crossings  at  ou  =  ^  1/2  and  ±  3/2.  A  sketch  of  the  desired  rational 

function  R  K,,  is  shown  in  Figure  9. 
e  N 

From  an  inspection  of  Figure  9,  and  Q_  can  be  easily  determined. 
Placing  these  quantities  in  (2.13)  gives  the  rational  function 


e  N 


N-2 

N/2  2  ^2 

(a2N+l)  -  (2i-l)^]  +  (a)^-4i^) 

i=l  i=l 


N/2 


N-2 

2 


;  N  even 


i=l 


(a^j^+l)  jQ  [U3^  -  (2i-l)^]  -  (a2jj-l)a)^(uj^-b^)  (u)^-4i^) 

i=l 


(3.10) 


N-1 

2 


N-1 

2 


2  *^2 
(a  +1)  (oJ^-b^)  n  [U)^  -  (2i-l)^]  +  (a2jj-l)uJ^  H  (tu^-^i^) 


e  N 


i=l 


i=l 


N-1 

2 


N-1 

2  2 
2  ..2. 


N  odd. 


(a^j^+l)  (cu  -h'  )  "  (2i-l)  ]  -  n  ^ 

i=l  i=l 


,  Letting  u)  =  1/2  and  setting  the  numerator  equal  to  zero  gives  the 

following  equation: 


^^2N  (1  -  4bf) _  rr  f  1  +  4i 

^^2N  ■  (2N  -  1)^  (2N  -  3)^  iii  \  ' 


N-2 

2 

n 


=  0  ;  N  even  (3.11) 
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Figure  9. 


A  rational  function  which  approximates  cos  tt  w 
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R  K 
e 


5 


-160(1) 


10 


+  6480(1)^  -  85290uj^  +  417460(1)^  - 
2u)^°  +  24600®  +  2740(1)^  +  26577u)^ 


63697500^  +  134460 
+  134460 


(3.14) 


e  D 


240ou^^ 


1  ^  Q  a  A  o 

-  16920(1)  +  42133500  -  4501140a)  +  1995637500  -  29216790(1)  +  6125625 

nO  ID  ft  A  A  Q 

2(1)  -  20(1)  +  47100  +  4780(1)  +  88727u)  +  1011690(1)  +  6125625 


A  check  of  the  accuracy  of  the  sinusoidal  approximation  at  o)  =  1/4 
gives  an  error  of  .248,  .0184  and  .00296  percent  for  N  =  2,  3  and  4  respectively. 

This  section  on  real  part  approximation  will  be  completed  after  a 

novel  technique  is  presented  for  obtaining  linear  phase  all-pole  transfer  functions. 

From  the  theory  of  the  Hilbert  transform,  if  the  real  part  of  a  network  function 

13 

is  cosinusoidal  then  the  imaginary  part  will  be  sinusoidal.  Hence,  the  phase 


0 


-1  im  part 

tan  - - - 

real  part 


tan 


sin  n  (1) 
cos  n  u) 


n  (D 


will  be  linear. 


(3.15) 


The  one  drawback  to  the  above  is  that  the  sine  is  the  Hilbert  trans¬ 
form  of  the  cosine  only  if  the  latter  is  defined  for  ail  frequencies.  If,  how¬ 
ever,  the  real  part  is  cosinusoidal  over  the  finite  frequency  interval  U)q) 

and  zero  elsewhere,  then  the  imaginary  part  will  be  roughly  sinusoidal  in  the 
vicinity  of  a  =  0.  The  reasoning  for  this  last  statement  is  best  explained  by 
noting  the  Hilbert  transform  of  the  cosine 

CO 

Im  part  =  +  —  f  ^  .  (3.16) 

n  J  §  -  U) 

—00 


This  integral  relationship  indicates  that  the  behavior  of  the  imaginary  part, 

at  a  frequency  u)  ,  is  due  primarily  to  the  integrand  within  the  vicinity  of 
°  1 

U)  .  The  modulation  of  the  cosine  function  by  ■= -  is  so  slow  outside  this 

o  5  -  u) 

vicinity  that  the  positive  and  negative  areas  nearly  nullify  each  other. 

The  real  part  function  developed  in  Example  (3.2)  is  near  cosinusoidal 

in  the  interval  (-N,  N)  and  approaches  zero  elsewhere.  Hence,  using  the  well 

30 

known  Brune-Gewertz  technique,  the  function  G  (s)  can  be  obtained  from 
2 

R  G„(a)  ) .  This  method  involves  the  factorization  of  the  denominator  polynomial 
e  N 
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and  the  solution  of  several  simultaneous  linear  equations;  these  tasks  make 
it  a  very  laborious  problem  for  high  order  real  parts.  To  overcome  these 
difficulties,  a  digital  computer  program  was  written  to  perform  the  necessary 
operations.  The  resulting  Gj^(s)  obtained  by  this  program  are  given  below  for 
N  =  2-6. 


%  ~  2 


-1.08719s  +  .707107 


s  +  1.18721s  +  .707107 


*^3  3 


1.87170s  -  2.16633s  +  1.73205 


s  +  2.04823s  +  3.09762s  +  1.73205 


""4  =  4 


3  2 

-2,64461s  -f  4.34677s  -  10.1194s  +  6.36396 


s'  +  2.89390s^  +  8.18732s^  +  10.4980s  +  6.36396 


(3.17) 


,  _  3.41295s^  -  7.25399s^  +  33.0873s^  -  39.2985s  +  28.4605 

^  s^  +  3,73461s^  +  16.9736s^  +  34,8719s^  +  47.6806s  +  28.4605 


^6  ~  6 


-4. 179258^  +  10.8879s^  -  816213s^  +  139.974s^  -  252.955s  +  159.099 


s"  +  4.57310s^  +  30.45663“*  +  86.7753s^  +  200.779s^  +  258.945s  +  159.099 


Though  it  isn't  immediately  apparent  without  factoring  the  numerator, 

all  the  zeros  of  G.,(s)  are  in  the  right  half  of  the  s-plane.  Furthermore,  the 
N 

phase  of  a  right-half  plane  zero  at  is  indistinguishable  from  the  phase 

of  a  left  half  plane  pole  at  -a^  +  Therefore,  if  an  all  pole  transfer 

function  G„(s)  is  desired,  it  can  be  obtained  from  G„(s)  =  N(s)/D(s)  by  replacing 
N  N 

s  with  -s  in.  the  numerator  N(s),  and  then  placing  N  in  the  denominator.  Hence, 

-  A 

G„(s)  =  — — ,  V  where  A  is  a  constant. 

N  N(-s)  D(s) 

From  the  so  obtained,  the  group  delay  response  was  calculated  and  plotted 

as  shown  in  Figure  10.  From  this  plot  it  is  concluded  that: 

a.  The  group  delay  approximates  rr  as  predicted  in  (3.15). 

b.  The  approximation  is  better  near  the  origin  as  predicted. 

c.  The  approximation  is  better  near  the  origin  for  higher  values  of  N. 
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d.  Tfib  group  delay  is  relatively  constant  within  the  6  db  passband. 

e.  The  number  of  group  delay  ripples  is  not  a  maximum. 

3.2  Imaginary  Part  Approximation 

The  imaginary  part  of  a  network  function  is  an  odd  rational  function 
of  the  real  variable  uj.  This  is  easily  shown  by  considering  the  odd  part 


Odd  H(s) 


^1^2  -  ""2^1 
2  2 
“2  -  ^2 


(3.18) 


of  the  network  function  of  (3.1).  It  is  now  clear  that  upon  letting  s  =  jo) 
the  above  will  be  an  odd  rational  function  in  the  real  variable  U)  multiplied 
by  j.  Hence,  it  is  the  imaginary  part. 

The  primary  property  of  the  imaginary  part  is  that  it  is  odd.  How¬ 
ever,  if  the  rational  function  of  Chapter  2  is  to  be  identified  as  an  imaginary 

30 

part,  it  must  also  possess  the  properties 

a.  the  numerator  polynomial  must  be  odd 

b.  the  poles  must  have  quadrantal  symmetry 

c.  the  degree  of  the  numerator  must  be  less  than  or  equal  to  the 
degree  of  the  denominator  plus  one,  equality  holding  if  H(s) 
has  a  pole  at  s  =  “. 

Several  specific  examples  will  now  be  given  wherein  an  odd  poly¬ 
nomial  is  approximated  by  an  odd  rational  function  to  give  an  equiripple  error. 

In  addition  it  will  be  seen  that  the  rational  function  so  obtained  also  satisfies 
the  above  listed  properties  of  an  imaginary  part. 

Example  (3.4) 

An  odd  rational  function  is  to  be  found  which 

a.  approximates  the  polynomial  =  u)  in  the  interval  (7IO,  10) 

b.  has  a  peak  error  of  unity  in  the  passband 

c.  has  a  4th  degree  denominator  and  1st  degree  numerator,  and 

d.  satisfies  the  properties  of  an  imaginary  part  of  a  network  function. 


Such  a  rational  function  is  illustrated  in  Figure  11. 
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Solution: 

Both  the  rational  function  and  the  polynomial  are  odd,  hence  the 
rational  function  can  be  determined  from  the  simplified  form  of  (2.20).  By 
inspection; 

2  2 

Odd  (u)  -  1)  (u)  +  10)  (u)  +  (JU-)  (u)  -  U)  ) 

H(a))  - - 2-^ - -  (3.19) 

Ev  (oo  +  10)  (w  +  (U  )  (u)  -  tu  ) 

A  X 


H(a))  = 


(3.25) 
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This  function  satisfies  the  properties  of  an  imaginary  part,  and  a 
Z(s)  could  be  found  with  the  above  H(u))  as  its  imaginary  part  if  desired. 

Example  (3,5) 

Whenever  high  order  polynomials  are  to  be  approximated  by  rational 

functions  whose  numerator  degree  n  is  less  than  the  degree  of  the  denominator 

m,  several  complicated,  non-linear,  simultaneous  equations  have  to  be  solved. 

This  example  is  chosen  to  illustrate  this  point.  In  the  interval  (f-10,  10)  an 

3 

odd  rational  function  G  is  to  be  obtained  such  that  it  approximates  =  o) 
in  an  equiripple  manner.  Two  solutions  are  obtained.  One  gives  a  maximum 
error  ripple  of  one  with  the  degrees  m  and  n  equal  to  4  and  3  respectively. 

The  othar  solution  considers  the  problem  with  m  and  n  equal  to  4  and  1  respec¬ 
tively.  The  anticipated  approximation  is  shown  in  Figure  12. 

Solution: 

By  inspection  of  Figure  12  and  using  (2.20) 

3  2  2  2 

Odd  ((JO  -€)  ((u  4  10)  (cu  +  a»~)  (o)  -  to.)  -  oUg) 

G((jd)  = - 5 - g - 2 - 

Ev  (cjo  +  10)  (ou  +  (JDg)  ((JO  -  (U^)  (UJ  -  CUg) 

Note  that  as  (3.26)  now  stands,  m  =  6  and  n  =  9.  It  appears  as  though  fopr 
coefficients  have  to  be  set  equal  to  zero  in  order  to  satisfy  the  first  sub¬ 
case.  Upon  closer  examination  of  (3.26)  however,  it  is  seen  that  the  co- 
6  9 

efficients  of  oo  and  ou  are  identical,  thus  reducing  the  number  of  simultaneous 
equations  to  three. 

The  polynomial  coefficients  of  the  numerator,  N((d)  and  denominator, 
D({ju),  are  obtained  in  a  straightforward,  but  lengthy,  manner.  They  are 

N(ou)  =((104-24>(U?' +)  (10a^+20b+2c+2ab-e)ou^  +  [l0b^+20ac+2bc-e  (20a+a^+ab)  ]uu^ 

+  [l0c^-€(20c+20ab+b?+2ac) ]w^-  C  (20bc+c^)u)  (3.27) 

6  2  4  2  2  2 

D((i))  =  (10+2a)ou  +  (10a  +20b+2c+2ab)(jD  +  (10b  +20ac+2bc)(JU  +  10c 
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It  is  noted  that  if  €  is  known  G((ji))  can  be  obtained  directly  from 
b,  and  c  without  having  a  knowledge  of  0)2,  and  o)^.  If,  however,  these 

values  are  desired,  the  cubic 

3  2 

B  +aB  +  bB  +  c  =  0  (3.30) 


can  be  factored,  the  zeros  of  which  are  then  a)_,  U)  ,  and  -  a)„. 

Letting  6=1,  the  equations  (3.29)  can  be  used  to  obtain  a,  b,  and 
c.  From  these  and  (3,27),  the  coefficients  of  the  U)  terms  can  be  calculated. 
The  G  so  obtained  is 


G(uj) 


152474.910)^  +  111640. 59u) 

-  175.60240)^  +  161305.3 


(3.31) 


The  frequencies  of  maximum  error  can  be  obtained  from  (3;30).  They  are 

0),  =  2.24838 
4 

0)^  =  6.26492  (3.32) 

=  9.01654 

This  is  one  of  the  required  solutions  to  the  approximation  problem. 
The  other  is  more  difficult  to  obtain  in  that  the  ripple  factor  £  must  be  used 
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to  set  the  coefficient  of  u)  equal  to  zero.  This  particular  coefficient  is 
a  function  of  a,  b,  c,  and  € .  Since,  however,  the  dependence  of  a,  b,  and 
c  upon  e  is  known  from  (3.29),  this  coefficient  can  be  determined  solely  as 
a  function  of  € .  With  this  coefficient  set  equal  to  zero,  and  after  the 
usual  algebra,  the  following  fourth  order  equation  in  £  results; 


(3.33) 


Factoring  the  above,  the  possible  €'s  are  found  to  be 


€  =  33.26862;  250;  1321.9980;  -355.26665  . 


From  these  roots  it  is  seen  that  even  the  smallest  value  of  ripple  is  com¬ 
paratively  large.  Using  the  smallest  value  the  corresponding  G  is 


G(U))  = 


_ 205383. 5u) _ 

(U^  -  198.02523u)^  +  11927.04 


(3.34) 


The  maximum  error  then  occurs  at  the  frequencies 


=  2.94838 
4 

gUg  =  7.25754 
U)^  =  9.30916  . 


(3.35) 


If  desired,  the  Brune-Gewertz  method  could  be  used  to  obtain  a 
function  Z(s)  which  has  either  of  the  above  G(a))  for  its  odd  part. 

3.3  Magnitude  Approximation 

The  properties  of  the  magnitude  function  are  easily  obtained.  From 
the  network  function  H(s)  defined  in  Section  3.1,  the  magnitude  squared 
function  is  found  to  be; 


22  “l^  ■  ^1^ 

|H(a)  )|  =  H(j(ju)  H(-jU))  =  - 5 - 5 

«2  -  ^2 


(3.36) 
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From  this  it  is  concluded  that  any  even  rational  function  of  uo  having  the 
properties  that^"^ 

a.  all  poles  and  zeros  have  quadrantal;  symmetry 

b.  any  pole  or  zero  on  the  real  frequency  axis  is  of  even  multiplicity 

c.  the  function  is  positive  for  all  frequencies 

can  be  used  for  a  magnitude  squared  function. 

The  magnitude  squared  function  most  often  used  for  filters  is  con- 
6 

veniently  expressed  as 


!h((u^) 


,  2  k  .  2. 

1  +  a  (u  F(cu  ) 


;  k  =  0,  1 


(3.37) 


k  2 

where  a  is  a  positive  real  constant  less  than  one  and  the  function  u)  F(oo  )  is 

chosen  in  such  a  manner  that  it  is  small  (usually  of  unity  peak  value)  in  the 

k  2 

passband  and  large  in  the  stopband.  This  cu-  F((ji)  )  will  then  give  a  magnitude 

14 

squared  function  which  approximates  the  ideal  square  magnitude  function. 

k  2 

In  Example  3.6  an  tu  F(a)  )  is  found  by  the  method  of  Chapter  2  such 

23 

that  the  ideal  low  pass  magnitude  function,  defined  as 

1  for  -1  <  0)  <  1 


0  elsewhere 


is  approximated  in  an  equiripple  manner. 

Example  3.6 

k  2  ^ 

Several  u)  F((i)  )  are  to  be  found  such  that  the  degenerate  polynomial 
=  1/2  is  approximated  with  a  ripple  factor  €  =  1/2  in  the  interval  (-1,  1). 
Outside  of  this  interval  the  approximating  function  is  to  be  large. 

Solution: 

To  satisfy  the  last  statement  of  the  problem  demands  that  n  >  m. 

Hence,  depending  upon  whether  the  equal  or  greater  than  sign  is  used,  the 

approximation  fits  either  Case  II  or  III.  For  the  present  discussion  assume 

the  maximum  number  of  ripples; then  a  desirable  rational  function  would  appear 

as  shown  in  Figure  13.  Using  the  notation  of  Figure  13,  the  rational  function, 
2 

call  it  f„(UJ  )  for  brevity,  is 
N 
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N-I 

1 

0)^  n 

i=l 


(3.38) 


N-1  ■  N-1 

2  2 
2  TT  /  2  2  ^2N''^  ,2  ,,  TT  ,  2 

^  2i-l>  -  il  -  “^21  ) 

.  ,  2N  ,  , 


N  odd 


i=l 


i=l 


N 

2 


1 — r  /  2  2 

n  2i-i)2 


-p  /•  2,  i=l 


N 

2 


N-2 

2 


n,  2  2  ^  ^2n"‘’  2  2  rr  p  2  2  ^ 

2i-l^  -  IT -  -  1)  il  (U)  -  CD  2.) 


i=l 


2N 


i=:l 


N  even 


where 


b  = 


0  if  m  <  n 


1  if  m  =  n 


To  obtain  the  classic  Chebyshev  polynomials  the  oj^'s  are  so  chosen 

that  all  the  poles  of  (3,38)  occur  at  cu  =  infinity.  This  means  that  b  =  0 

2i 

and  the  N-1  denominator  coefficients  of  id  ,  i=l,...,  N-1,  are  equated  to 
zero.  This  gives  N-1  equations  and  as  many  unknowns.  Hence  all  ud^'s  are 
uniquely  specified  by  the  process.  An  illustration  of  the  method  for  N=4 
is  given  as  follows;  Expanding  the  denominator,  the  polynomial 
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(D  (1  +  20)2  ”  ^^1  ~  ^^3  ^  ^  “^3  ‘^'^1  ^^3  "  ^^^2  “  *^2  ^  (3.39) 

2,2„42„4  2^  44 

+  0)  (cu^  “  2u)^  tUg  -  20)^  )  +  OJ^  CDg 


results.  The  3  simultaneous  equations  are  then 


2  2  2 
1  +  20)^  =  2(u)^  +  u)^  ) 


2  4  2  2  2  2 

2w^  +  U)^  =  (cUj^  ‘"s  ^ 


4  „  2  2  ,  2  2 

“2  =  ^1  ‘"3  ^‘"l  ^  "3 


(3.40) 
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Letting  x  =  o)^  +  ID^  ;  y  =  the  solutions  x 

obtained.  Using  these  results 


1,  and  y  =  1/8  are 


o  4  2  2 

i'^(u)  )  =  (8u)  -  8u)  +1) 


g 

which  is  recognized  as  the  square  of  the  Chebyshev  polynomial  of  4th  degree. 

Another  classical  design  is  the  elliptic  filter  with  an  equiripple 

6 

magnitude  in  both  the  pass  and  the  stopbands,.  It  can  be  shown  that  for  such 
2 

a  filter  the  )  function  has  poles  and  zeros  which  are  reciprocally 

related^  i.e.,  if  oi^  is  the  proportionality  constant, 

U)  .(ju  .  =  a>  (3.41) 

pi  zi  c 

21  2 

where  O)  .  (o)  .)  is  the  location  of  the  i^h  pole  (zero).  A  sketch  of  f„(u)  ) 
pi  zi  N 

is  shown  in  Figure  14. 


The  elliptic  functions  for  N=2  and  3  can  be  obtained  rather  easily 
from  (3.38).  Letting  b  =  1  and  0  for  these  two  respective  cases,  the  required 
polynomial  equivalences  are 


2  2. 

(OD  -  )  - 


a^-lN 


2/2  ..  1 
OU  (OD  -  1)  =  -7 


'  2  '^c 

(U - ; 


N  =  2 


(3.42) 


22  2^  2  2  2^1^2 
U)  (U)  -  )  -  (uj  -  1)  (UD  -  U)2  ^  ~  17^  ['^ - 1 

mV  u), 


;  N  =  3  . 


Equating  coefficients  and  solving  for  the  unknowns  in  terms  of  (ju  and  M 

c 

respectively,  gives  the  following  solutions 


2  2 
2  (u)  -  ) 


(3.43) 


M  = 


1) 


2 


where 
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l/ou  . 
c 


Similarly, 


2,  2  2^ 

„  U)  (OJ  -  (JU  ) 

=  r - ^ 

2 

2  *^0  2 

(CU - 2) 

‘"l 


(3.44) 


where 


2 


3M  -  1 
M  +  1 


2 


U) 


2 


2  4 

OU  =  MD- 
c  2 

=  f  (UJ  . 
3  C 


Both  solutions  agree  with  those  given  by  Papoulis  ,  who  arrived  at 
them  in  a  somewhat  different  manner. 

For  higher  orders  of  N,  extensions  of  the  above  procedures  would  lead 

to  several  non-linear  simultaneous  equations;  however,  because  this  problem 
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has  been  given  a  more  elegant  solution  in  the  literature  '  it  will  not  be 
pursued  further , 


3 . 4  Constant  Group  Delay  Approximation 

The  group  delay,  t,  of  the  all-pole  transfer  function 


H(s)  = 


“2  ^  ''2 


N 

n 

i=l 


R  s.  <  0 
e  1 


(s  -  s^) 


(3.45) 


where  M„  +  N  are  the  even  and  odd  parts  respectively  of  a  Hurwitz  polynomial, 

^  ^  14 

can  be  expressed  in  one  of  two  ways.  The  first 


T 


Where  '  Indicates 
differentiation 


(3,46) 


1 

J 


^2^2’ 


s=j(u 


is  the  polynomial  representation. 


The  second 


29 


1 


s  +  s . 
1 


] 


S=j(JU 


I 


s=ju) 


is  in  terms  of  the  left-half  s-plane  pole  locations  of  H(s). 

Using  analytical  continuation  such  that  T  is  valid  for  all  complex  s, 

29 

the  following  properties  of  t  can  be  deduced; 

a.  t(s)  has  two  zeros  at  s  =  oo 

b.  t(s)  is  an  even  function 

c.  t(s)  is  positive,  real  and  bounded  for  s  =  juj. 

The  necessary  and  sufficient  conditions  that  an  even  rational  function  of  s  be 
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a  group  delay  function  of  an  all  pole  transfer  function  are: 

a.  The  poles  have  quadrantral  symmetry,  and  there  are  no  poles  on 
the  ju)  axis, 

b.  The  residues  of  left-half  plane  poles  are  1/2,  right  half  plane 

2 

poles  -1/2.  Equivalently  if  s  is  replaced  by  p,  the  resiude  of 
any  pole  is  the  square  root  of  the  pole  location. 

From  the  foregoing,  it  is  now  possible  to  pick  one  of  the  equiripple 
rational  functions  developed  in  Chapter  2  and  force  it  into  being  a  group  delay 
function.  At  the  outset,  it  is  noted  that  the  rational  function  chosen  must 
(1)  be  an  even  function  and  (2)  have  a  numerator  degree  of  two  less  than  the 
denominator.  Condition  (2)  and  the  fact  that  a  constant  is  being  approximated 
fix  the  approximation  to  Case  IV.  The  degree  of  the  denominator  will  be  2N, 
hence,  from  Table  1  the  total  number  of  arbitrary  constants  is  n  +  4  =  2N  +  2. 
Two  of  these  constants  are  specified  by  fixing  the  polynomial  =  c^  and  the 
relative  ripple  €.  Also  from  Table  1,  the  maximum  number  of  peak  deviations  is 
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2N  +  1.  Because  the  function  is  even,  one  extreme  point  will  always  occur  at 
the  origin;  the  others  will  be  symmetrical  about  u)  =  0.  Thus,  the  total 
maximum  number  of  arbitrary  peak  deviations  is  N,  but  this  is  exactly  the  num¬ 
ber  of  arbitrary  constants  yet  to  be  specified.  Therefore,  it  can  be  con¬ 
cluded  that 

a.  The  polynomial  =  c^  can  be  made  any  convenient  value  (for  this 
problem  let  c^  =  1) . 

b.  The  ripple  factor  e  is  arbitrary. 

c.  The  N  pole  locations  can  be  fixed  by  varying  the  N  unspecified 
peak  deviation  frequencies  until  the  residue  condition  is  met. 

The  resulting  left  half  plane  poles  of  the  equiripple  rational 
function  can  then  be  used  to  obtain  a  transfer  function  H(s) 
with  an  equiripple  group  delay. 

d.  The  equal  ripple  group  delay  will  look  similar  to  the  illustration 
of  Figure  15. 

Using  the  notation  of  Figure  15  and  the  rational  function  in  (2.14b), 
the  proposed  group  delay,  as  a  function  of  ou  is: 


N-1 

2 


n 


T(u))  =  1  +  e  -  2e 


i=l 


N-l 

2 


N-l 

2 


N  Odd 


n,2  .2  l-e,2  2.  77,2  2  , 

-  ‘“2i>  -  m  -  “'i  >  ii  -  “^Bi  i:  1^ 


i=l 


i=l 


(3.48) 


N 

2 


T(u))  =  1  +  e  -  2e: 


n -  ‘“2i'> 

i=l 


N-2 

2,^  1  -  e  2  ,  ,2  2,  .  2 


•;  N  even. 


—  -z  a.  i  -  t  a  .  ^  y-r  .  ^ 

n  ■  ‘“zi  >  -  ITT  (u)  -  )  n  (u)  -  0^2,  ^ 


2  2 


1) 


i=l 


i=l 
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2  2 

Using  analytic  continuation,  i.e.,  letting  s  =  -  lu  ,  and  then,  since  T  is  an 

2 

even  function,  letting  s  =  p  the  above  becomes 


2 

where  u.  =  0U.  ;i=l,  2,  ....  N. 

1  1  ’  ’  ’  ’ 

It  is  interesting  to  note  at  this  point  that  the  results  of  Ulbrich 
29 

and  Piloty  could  be  duplicated  by  letting  u^^  =  1  and  letting  the  average 
delay  be  an  unknown  quantity.  The  ripple  factor  £,  under  these  circumstances, 
would  be  an  absolute  delay.  Because  of  the  transformation  involved,  Ulbrich 
and  Piloty  could  not  be  this  flexible. 

2 

Abele  did  not  use  an  explicit  rational  function  in  his  work  but 
imposed  constraints  on  the  numerator  and  denominator  polynomials  of  (3.46). 

These  constraints  appear  in  a  form  similar  to  the  numerator  and  denominator  of 
(3.48). 

Prior  to  the  publication  of  Abele 's  work  this  author  worked  out  closed 

form  solutions  to  the  equiripple  group  delay  for  n  =  1  and  2,  using  an  approach 

16 

similar  to  that  of  Abele.  Helman  ,  using  Bernstein's  Chebyshev  rational 
function,  also  obtained  solutions  for  n  =  1  and  2f.  In  Helman 's  technique,  how¬ 
ever,  it  is  necessary  to  find  the  aeros  of  a  high  order  polynomial  to  obtain 
the  desired  solution. 
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where 
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Equating  the  coefficients  of  (3.51)  to  those  in  (3.52)  and  solving 
the  resulting  simultaneous  equations,  the  solutions  lor  n  =  1  and  2  are  found 
to  be 

1  / - 

COi  =-^y2£(l  -  £) 

(3.53) 

1 

^o  “  1  +  € 


n  =  2; 


Obviously,  the  equations  for  n  greater  than  two  become  sO’  complicated 
that  a  closed  form  solution  is  no  longer  possible.  Numerical  answers,  on  the 
other  hand,  can  still  be  obtained  by  using  a  high  speed  digital  computer  to 
perform  the  many  operations  required  in  an  interative  solution. 

There  are  many  approaches  which  can  be  used  in  solving  non-linear 
1  17 

simultaneous  equations.  ’  One  method  is  the  Newton-Raphson  technique.  This 
procedure  has  an  advantage  over  many  others  in  that,  if  the  initial  guess  is 
within  the  region  of  convergence,  the  convergence  is  extremely  rapid.  A  dis¬ 
advantage  of  the  method  lies  in  the  fact  that  the  region  of  convergence  lor 
many  problems  of  this  type  is  small. ^  For  the  problem  at  hand  the  u)^'s  will 
be  iterated.  This  means,  if  the  Newton-Raphson  procedure  is  used,  that  the 
initial  guesses  must  be  known  reasonably  accurately.  Such  an  estimate  can  be 
made  very  easily  if  it  is  reasoned  that  the  uo^'s  lor  increasing  N  will  follow 
a  definite  pattern.  Furthermore,  this  pattern  has  already  begun  with  a  closed 
form  solution  for  n  =  1,  2. 


57 


Briefly,  the  method  of  solution  undertaken  is  as  follows: 

a.  From  a  knowledge  of  lower  order  solution,  a  set  of  (o^'s  is 
picked  to  correspond  to  a  large  ripple  factor  €  (it  is  found 
experimentally  that  the  region  of  convergence  is  larger  lor 
large  €) . 

b.  Using  these  values  of  (U^,  the  denominator  of  (3.49)  is  factored 
and  the  residues  obtained. 

c.  Knowing  that  the  residues  squared  less  the  root  locations  must 
be  zero  gives  N  equations,  F^,  to  be  satisfied  in  terms  of  the 
N  unknowns  U).  ,  i=l,  ...  N. 

d.  The  partial  derivatives  3F^/3od^  are  obtained  numerically  in  the 
approximate  form 

3F.  F,(u).  +  6u).)  -  F,(a).) 

_ J  _  J  1  1  J  1 

where  6  is  a  small  constant. 

e.  The  N  simultaneous  equations  (Newton-Raphson) 

hV . 

A(ju.  =  -F.  ;  j  =  1,  ...  ,  N  (3.55) 

dU).  1  J 

i=l 

are  solved  for  the  unknowns  Au).  . 

1 

f.  The  initial  values  of  ou.  are  replaced  by  u).  +  Au).  . 

1  i  1 

g.  The  above  is  repeated  (usually  from  4  to  8  times)  until  the 
(ju^'s  are  made  as  small  as  the  numerical  accuracy  of  the  com¬ 
puter  allows.  The  error  commonly  obtained  is  from  one  part 

10  7 

in  10  to  one  part  in  10  for  N  ranging  from  3  to  10. 

h.  When  the  set  of  solutions  ou.  is  known,  it  is  then  used  for 
the  initial  guesse,  and  the  ripple  factor  €  is  given  a  smaller 
value. 

The  pole  locations  of  the  all  pole  equiripple  group  delay  filter 
resulting  from  the  above  program  are  tabulated  in  Table  2  of  the  appendix  for 
N  =  3,  ...  ,  10,  each  with  €  =  0.005,  0.01,  and  0.05.  For  comparative  purposes. 


TABUS  2 


Root  locations  of  an  N  pole  (N  =  2  to  10)  equlripple  group  delay 
filter  for  3  values  of  ripple,  €.  The  normalization  is  such  that  the  mean 
low  frequency  delay  is  unity. 
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3 


4 


5 


6 


7 


8 


9 


10 


e  = 

o;oo5 

e  = 

0.01 

e  = 

0,05 

real  part. 

imag.  part 

real  part 

imag.  part 

real  part 

imag.  part 

-1.4021451 

i. 85237134' 

-1.3628732 

*.94649670 

-1.2039728 

*1.0416954 

-1.9088240 

-1.5693900 

0.0 

•*^2.0151812 

-1.8061331 

-1.5022964 

0.0 

*2.0802275 

-1,4703023 

-1,2765428 

0.0 

*2,2963311 

-2.1050989 

-1.6498925 

^1. 1051383 
*•3.2352073 

-1.9643081 

-1.5669472 

*1.1459939 

*3.3388487 

-1.5460122 

-1.3079526 

*1.2621897 

*3.6530753 

-2.3081926 

-2.1996486 

-1.6951436 

0.0 

*■2.3365664 

*4.5363797 

-2.1262469 

-2.0381080 

-1.6028093 

0.0 

*2.4144421 

*4.6721443 

-1.6218304 

-1.5791674 

-1,3250626 

0.0 

*2.6242135 

*5.0642095 

-2,4057414 

-2.2532646 

-1.7235697 

*1.2339268 

*3.6461522 

*5.8905328 

-2.2016824 

-2.0794685 

-1.6252314 

*1.2704024 

*3.7547734 

*6.0529615 

-1.6550387 

-1.5974133 

-1.3357252 

*1.3628473 

*4.0389428 

*6.5084283 

-2.5044409 

-2.4611483 

-2.2872444 

-1.7428862 

0.0 

*2.5460921 

*5.0065049 

*7.2812168 

-2.2777117 

-2.2440387 

-2.1055654 

-1.6404479 

0.0 

*2.6127150 

*5.1407019 

*7.4661369 

-1.6882400 

-1.6733428 

-1.6088716 

-1.3429725 

0.0 

*2,7783913 

*5.4857180 

*7.9747231 

-2.5604843 

-2.4963531 

-2.3105097 

-1.7567883 

*1.3128118 

*3.9086731 

*6.4018708 

*8.6983493 

-2.3203970 

-2.2708293 

-2.1234089 

-1.6513970 

*1.3427945 

*4.0003415 

*6.5577090 

*8.9025775 

-1.7065422 

-1.6848570 

-1.6167055 

-1.3482052 

*1.4157328 

*4.2258593 

*6.953096 

*9.4567365 

-2.6168521 
-2.5961134 
-2.5205210 
- .23273558 
-1.7672335 

0.0 

*2.6761673 

*5.3058951 

*7.8226335 

10.135398 

-2.3632467 

-2.3474107 

-2.2891886 

-2.1363241 

-1.6596318 

0.0 

*2.7309954 

*5.4187559 

*7.9971129 

*10.356477 

-1,7248403 

-1.7180599 

-1.6927410 

-1.6223879 

-1.3521532 

0.0 

*2.8634253 

*5.6946244 

*8.4373817 

*10.950485 

-2.6526863 

-2.6205920 

-2.5380595 

-2.3400882 

-1.7753412 

*1.3635814 

*4.0741370 

*6.7281919 

*9.2625868 

*11.587925 

-2.3903646 

-2.3659368 

-2.3025079 

-2.1460838 

-1,6660347 

*  1.3883671 
*4.1499642 
*6,8593217 
*9.4533690 
11.823921 

-1.7363559  *1.4477569 
-1.7259501  *4.3324074 
-1.6984663  *7.1785723 
-1.6266905  *9.9322817 
-1.3552374  *12.453322 
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Table  3  gives  (a)  the  zero  frequency  delay  (for  a  unity  6  db  bandwidth),  (b) 
the  ratio  of  the  3  db  to  6  db  attenuation  points,  (c)  the  ratio  of  the  equiripple 
approximation  interval  to  the  6  db  attenuation  points,  and  (d)  the  percentage 
overshoot  of  the  impulse  response.  Assuming  the  maximum  of  the  equiripple  group 
delay  transfer  function  occurs  at  zero  frequency  (this  assumption  is  true  for 
N  odd;  and  for  N  even  if  t  is  smaller  than  about  0.05),  the  numerator  constant 
of  the  transfer  function  can  be  made  equal  to  the  denominator  constant  and  the 
magnitude  characteristics  vs.  angular  frequency  can  be  computed.  A  plot  of  this 
characteristic  is  shown  in  Figure  16  for  N  =  3,  ...,  10  and  t  =  0.005.  The  cor¬ 
responding  attenuation  curves  for  0.005  <  0.05  are  but  slightly  different.  As 
t  becomes  larger  (s^O.l)  the  magnitude  response  departs  from  that  given  in  Figure 
16.  Two  major  differences  are  noted:  (1)  the  individual  poles  in  the  transfer 
function  become  apparent  by  the  lumpy  nature  of  response,  and  (2)  the  high  fre¬ 
quency  attenuation  is  a  few  db  greater. 

For  small  €,  the  attenuation  curves  are  almost  identical  to  those  of 

28 

the  corresponding  Bessel  or  maximally  flat  delay  filter  up  to  about  the  10  db 
attenuation  points.  For  higher  frequencies  the  attenuation  becomes  larger. 

The  applications  of  constant  group  delay  filters  are,  primarily,  for 
the  transmission  of  pulsed  information.  Thus,  the  impulse  response  of  such  a 
filter  is  of  prime  interest.  Two  families  of  curves  are  given  in  Figures  17 
and  18  for  the  impulse  responses  (N  =  3,  ...,  10)  with  t  =  0.005  and  0.05 
respectively.  From  these  responses  the  following  conclusions  can  be  drawn: 

a.  The  pulses  are  nearly  symmetrical  about  the  time  t  =  T^.  This  is 
particularly  true  for  large  N. 

b.  The  impulse  response  maximum  is  almost  independent  of  N. 

c.  The  overshoot  decreases  with  an  increase  in  N. 

d.  The  curves  are  not  plotted  to  sufficient  accuracy  to  observe 
the  peak  of  the  over  or  under  shoots.  However,  from  the  com¬ 
puter  data  it  is  noted  that  for  N  >  7  the  over  or  undershoot  ap¬ 
proaches  a  constant  of  0.2  and  2.0  percent  for  €  =  0.005  and 

€  =  0.05  respectively.  The  figures  predicted  in  the  introduction 
were  about  0.1  and  1.0  percent  respectively. 

The  reason  the  overshoot  and  asymmetry  is  greater  for  small  N  is  that 
the  group  delay  approximation  interval  is  not  as  large  as  the  passband  of  the 


TABLE  3 


With  a  unity  6  db  radian  frequency  bandwidth  (a)  the  low  frequency 
mean  delay,  (b)  the  3  db  bandwidth,  (c)  the  equiripple  group  delay  approxi¬ 
mation  interval,  and  (d)  the  impulse  response  overshoot  (€  =  0;005,  €  =  0.05) 
are  given. 

€  =  0.005 


N 

T 

O 

‘“3db/“6db 

‘“T/“^6db 

8^ 

2 

1.97316 

.701656 

.389385 

3 

2.51335 

.742110 

.  .715918 

2.3 

4 

3.11611 

.731226 

.962641 

2,9 

5 

3.62909 

.721310 

1.18577 

2.5 

6 

4.07041 

.718634 

1.38658 

1.7 

7 

4.45741 

.713878 

1.57773 

.9 

8 

4.82706 

.717055 

1.75023 

.4 

9 

5.16313 

.709500 

1.91449 

.2 

10 

5.48223 

.716544 

2,06791 

.2 

€  =  0.01 


2 

1.97566 

.707682 

.460990 

3 

2,53967 

.744943 

.793434 

4 

3.16442 

.727498 

1.03^60 

5 

3.64538 

.722188 

1.26177 

6 

4.11155 

.717549 

1.45456 

7 

4.47961 

.713319 

1.65060 

8 

4,86091 

,718993 

1.81672 

9 

5.18980 

.706297 

1,98182 

10 

5.50975 

.720700 

2,13315 

e  =  0.05 

2 

2,01019 

.737244 

.668501 

3 

2.63387 

.742351 

.992804 

3.1 

4 

3,33354 

.711244 

1.19497 

3.1 

5 

3.65745 

.736920 

1.47683 

2.2 

6 

4.29086 

.704986 

1.59647 

2.1 

7 

4.45445 

.722572 

1.86774 

2,1 

8 

4,99602 

.733237 

1.96240 

2.2 

9 

5.22731 

.691956 

2.16169 

2.2 

10 

5.59256 

.755409 

2.28950 

2.2: 

Figure  16.  Attenuation  of  the  Equiripple  Group  Delay  Filter  with  a  1/2  percent  ripple.  T 
izatlon  is  such  that  the  6  db  bandwidth  is  unity  (the  ideal  gausslan  filter  is 
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Figure  18,  Impulse  response  of  an  equiripple  group  delay  filtar  with  a  6  db  unity  bandwidth  and 
percent  ripple 
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filter.  Thus^  the  higher  frequency  components  arrive  out  of  phase  causing 
irregularities  in  the  impulse  response. 

3,5  Narrow  Band  Approximation  of  Linearly  Decreasing  Group  Delay 

This  chapter  will  be  concluded  after  the  formulation  of  the  equi- 

ripple^  linearly  decreasing,  narrow  band  group  delay  problem  has  been  invest!- 

20 

gated.  A  graphical  solution  to  this  problem  has  appeared  in  the  literature  , 
but  as  of  yet  an  analytical  one  has  not. 

Consider  an  all  pass  transfer  function 
^  (s  +  s. )  (s  +  s.*) 

t(s)  =  ]~f  - - — - - where  Res.  <  0  (3.56) 

11  (s  -  s.)  (s  -  s.*)  1 

i=l  ^  ^ 


with  an  even  number  of  complex  poles.  The  numerator  will  give  the  same  delay 
as  the  denominator.  Therefore,  the  total  delay  can  be  obtained  by  inserting 
a  multiplying  factor  of  two  into  (3,47),  giving 


T(jcJU) 


1 


s  +  s . 
1 


1 

s  +  s .  * 
1 


s=ja) 


(3,57) 


Letting 


T(ja)) 


=  1 


i=l 


^  „  (u3  +  00^)  +  j((u  -  (u^) 


_ 1 _ 

-CT^  +  j(U)  +  (h^) 


(3.58) 

_ 1 _ ■ 

-a^  +  j(cju  -  (Ju^)," 


The  narrow  band  approximation  is  obtained  by  noting  that  if  >.>  and  »  cu^ 

...  Rs  then  the  first  and  third  term  of  (3.58)  are  negligible  for  the 
N 

frequency  band  of  interest.  With  the  definitions 


0)^  =  a  constant  frequency  displacement 


W  =  UO  -  U) 

o 


U,  =  (1).  -  U) 

1  1  o 


V 


i 


a 

i 


u  +  jv 
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In  order  to  make  the  numerator  2  degrees  less  than  the  denominator 
requires  3  simultaneous,  non-linear  equations.  This  in  addition  to  the  residue 

condition  will  fix  the  N  +  3  unknown  constants  in  (3.60).  A  solution  of  this 
problem  has  not  been  attempted.  Because  of  the  large  number  of  equations,  it 
will  indeed  be  a  difficult  one. 


4.  A  numerical  procedure  fcr  obtaining  constant  group 

DELAY  WITH  EQUIRIPPLE  ERROR  AND  EXAMPLES 
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There  are  many  problems  of  an  equiripple  or  pseudo  equiripple  nature 
which  are  difficult,  if  not  nearly  impossible,  to  formulate  by  the  equiripple 
rational  function  of  Chapter  2.  Two  such  problems  are  considered  in  this  chapter. 
The  first  is  a  logical  extension  of  the  constant  group  delay  approximation 
undertaken  in  Section  3.4;  the  ripples  are  weighted  such  that  they  increase  with 
frequency.  The  second  is  the  determination  of  the  pole  locations  of  an  all¬ 
pass  network  to  compensate  the  non-constant  group  delay  characteristics  of  a 
sharp  cutoff  filter. 

The  proposed  method  of  solution  for  these  two  proDiems  utilizes  the 
rapid  convergence  properties  of  the  Newton-Raphson  iteration  technique.  The 
steps  taken  to  obtain  an  equiripple  or  pseudo  equiripple  solution  are  outlined 
below; 

a.  A  crude  solution  is  found  in  which  the  desired  number  of  ripples 
are  present.  This  solution  depends  upon  the  nature  of  the  prob¬ 
lem;  in  some  cases  it  is  easily  obtained,  in  others  it  is  more 
difficult . 

b.  If  there  are  N  arbitrary  unknown  parameters  available,  then  N 
equations,  F^,  i=l,...,  N  are  formed  by  setting  the  N  distances 
from  the  maxima  or  minima  to  a  known  reference  equal  to  zero. 

c.  The  maxima  and  minima,  if  not  known  analytically,  are  obtained 

17 

numerically  by  Newton's  method. 

d.  If  the  reference  mentioned  in  (b)  is  chosen  appropriately,  then 
after  one  iteration,  using  Newton-Raphson,  the  N  equations  are 
solved  with  sufficient  accuracy  that  the  reference  can  be  shifted 
towards  the  desired  ultimate  solution.  The  above  process  is  then 
repeated  until  the  final  reference  value  has  been  reached. 

The  two  problems  mentioned  above  will  now  be  used  as  examples  to 
illustrate  this  procedure. 

Example  (4.1)  : 

The  group  delay  can  be  made  to  approximate  a  constant  over  a  larger 
interval  by  increasing  the  ripple  factor.  This  Increase,  however,  causes  more 
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over  and  undershoots  to  be  present  (introduction) .  As  a  compromise  it  seems 
reasonable  that  the  group  delay  should  approximate  a  constant  best  where  most 
of  the  spectral  energy  falls,  or  where  the  magnitude  function  is  largest.  With 
this  type  of  weighted  error  the  approximation  interval  would  be  greater,  and 
there  would  be  less  distortion  in  the  impulse  response. 

The  problem  is  thus  two  fold,  (1)  to  pick  a  desirable  weighting 
function,  and  (2)  to  determine  the  pole  locations  of  an  all  pole  transfer 
function  such  that  the  peaks  and  valleys  of  the  group  delay  fall  on  the  boundary 
of  this  weighted  error  region. 

One  attractive  weighting  function  is 

(in  (ou/oi^)^  (4,1) 

error  =  ±  € 

where 

=  the  zero  frequency  error 

€  =  the  error  at  the  end  of  the  approximation  interval  , 

A  1 


With  this  weighting  function,  a  typical  group  delay  (n  =  7)  will  be  as  shown 
in  Figure  20,  Obviously,  a  crude  approximation  to  this  problem  is  already 
known,  i,e,,  the  equiripple  approximation  developed  in  Section  3.4,  The  N 
pole  locations  can  now  be  specified  in  terms  of  the  N  peak  deviations  occuring 
at  (ju  ,  ,,,,  uj^  by  means  of  the  N  non-linear  equations 

^  WtI. 


I'd  =T(^  .^^) 


i+1 

1  +  (-1) 


vs> 


=  0 
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and  the  maxima  and  minima,  the  latter  being  zeros  of 


d  T 


(U)^) 

dU)^ 


1 


s 

i 


+  s. 
1 


=  0  . 


(4.4) 


An  initial  choice  of  the  pole  locations  and  the  references  6^  and 

can  be  made  such  that  the  set  is  solved  to  a  reasonable  degree  of  accuracy 

in  one  iteration.  A  new  £  and  £  are  now  chosen  and  the  process  repeats 

X  ^ 

until  the  desired  solution  results.  Many  factors  dictate  the  number  of 
iterations  required  to  drift  the  poles  into  the  neighborhood  of  the  solution. 
Among  these  are;  (a)  the  required  overall  change  in  £j^  and  from  the 
starting  values,  and  (b)  the  pessimistic,  or  optimistic,  choice  of  the  incre¬ 
mental  changes  in  £^  and  A  pessimistic  choice  requires  many  more  itera-^ 

tions  than  necessary  and  an  overly  optimistic  choice  results  in  a  solution 
outside  the  region  of  convergence,  requiring  a  new  start.  Once  £,  and  £ 
reach  their  final  values,  continued  iteration  will  give  the  necessary  exacti¬ 
tude  to  the  pole  locations. 

A  computer  program  was  written  to  perform  the  above  operations.  Very 
few  iterations  were  required  for  the  final  answers  to  be  accurate  to  at  least 
8  significant  figures. 

A  typical  example,  with  N  =  10,  £^  =  0.01,  and  =  0.16,  has  the 
pole  distribution  shown  in  Figure  21.  For  comparison,  the  initial  pole  dis¬ 
tributions  of  the  equiripple  group  delay  filter  with  £  =  0.01,  and  a  Bessel 
filter  are  given  also. 

The  effect  of  the  weighted  ripple  approximation  on  the  impulse 
response  is  best  illustrated  by  a  low  order  example.  With  N  =  3,  £^  =  0.005 
and  =  0,08  the  impulse  response  is  computed  and  compared  to  that  of  the 
equiripple  group  delay  (€  =  0.005)  as  illustrated  in  Figure  22.  Two  differences 
are  noted  in  the  figure.  One  is  that  the  symmetry  is  improved  about  the  im¬ 
pulse's  maximum  due  to  the  increased  group  delay  approximation  interval.  The 
other,  increased  overshoot,  is  a  decided  disadvantage.  Both  of  these  differences, 
however,  are  small,  and  what  is  more,  these  differences  become  less  noticeable 
for  higher  N. 


Figure  22.  Comparison  of  the  equiripple  (solid)  and  the  weighted  ripple  impulse  response.  The 
6  db  bandwidth  is  unity  and  e  =  e  =  0.005,  €„  =  0.08 
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Example  (4.2): 

There  are  many  existing  procedures  for  compensating  the  non-constant 

V  X  2  15  27 

delay  characteristics  of  sharp-cutoff  filters.  >  >  >  Foremost  among  these 

methods  is  Darlington's  Chebyshev  polynomial  series  approach^^  which  many 

authors  have  used. This  method  requires  that  the  group  delay  be 

expanded  in  a  Fourier  series,  and  various  dominant  terms  in  the  series  are 

made  negligible  by  an  appropriate  choice  of  all-pass  network  constants. 

This  method  has  one  disadvantage,  i.e.,  the  number  of  Fourier  series  terms 

contributing  to  the  group  delay  error  is  unnecessarily  large  when  the  group 

delay  of  a  sharp  cutoff  filter  is  being  expanded  over  an  interval  that  extends 

12 

beyond  the  filter's  cutoff. 

This  example  is  used  to  illustrate  a  method  wherein  the  poles  of  an 
all  pass  compensating  network  can  be  found  to  give  an  overall  equiripple  group 
delay  characteristic  to  the  low  and  all-pass  filter  combination.  Further¬ 
more,  the  delay  can  be  made  equiripple  as  far  beyond  cutoff  as  is  desired. 

Before  the  problem  proper  is  started,  one  property  of  a  minimum 

phase  sharp  cutoff  filter  needs  to  be  investigated.  Guillemin,  using  the 

Hilbert  transform,  has  shown  that  the  group  delay  of  an  ideal  filter  (as 

defined  in  Section  3,3)  increases  in  a  monotonic  manner  until  it  reaches  a 

peak  at  or  near  the  filter's  cutoff.  Beyond  this  point  the  group  delay  de- 

14 

creases  monotonically  to  zero.  This  type  of  behavior  is  indeed  observed  in 

24 

the  classical  Butterworth  and  the  low  ripple  Chebyshev  filter  (for  larger 
ripples  the  group  delay  increases  but  in  a  non-monotonic  manner) .  The  group 
delay  of  a  typical  sharp  cutoff  lowpass  filter  is  shown  in  Figure  23.  In 
order  that  the  overall  filter  -  all-pass  combination  possess  a  constant  group 
delay,  the  shaded  area  must  be  added  by  the  all-pass  (an  all-pass  giving  neg¬ 
ative  group  delay  is  unrealizable).  An  equiripple  error  solution  to  this 
problem  can  now  be  obtained  in  the  following  sequence; 

a.  It  is  noted  that  the  narrow  peak  occur ing  at  the  cutoff 

is  due  primarily  to  the  pole  in  the  low  pass  filter  which  is 
nearest  the  cutoff  frequency,  and  also  nearest  the  jou-axis  of 
the  complex  s-plane. 

b.  An  equally  spaced  pole  distribution  parallel  to  the  jou-axis 
will  give  an  overall  group  delay  characteristic  that  is  pseudo 
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equiripple.  The  closer  the  poles  are  to  the  jou-axis,  the 

n  •  4-u  .  ,  14 

larger  is  the  ripple. 

c.  If  the  pole  distribution  of  (b)  is  chosen  such  that  it  includes 
the  pole  mentioned  in  (a),  then  an  "almost"  equiripple  group 
delay  can  be  achieved.  The  pole  locations  and  the  delay  char¬ 
acteristics  for  this  combination  are  shown  in  Figure  24. 

d.  In  the  work  that  follows  assume  that  an  even  number,  N,  of 
complex  poles  have  been  added  to  the  fixed  complex  pole. 

The  N  equations 


F.  =  T((U.)  -  t(uj.^2)  =  0 
i  =  1,  ...  N 


(4.5) 


will  be  of  an  even  number  and  can  be  solved  by  Newton- 
Raphson  to  determine  the  exact  locations  of  the  added  N 
poles  to  give  an  exact  equiripple  group  delay.  Once  N  and 
the  fixed  pole  location  have  been  specified,  no  further  control 
exists  over  the  ripple  €  or  the  average  delay  T^. 

e.  The  group  delay  for  the  above  equiripple  configuration  can  be 
written  in  the  form 


tCo)^) 


II 


5 

i 


2  2 
cu  +  s . 

1 


/  2, 

+  ) 


(4.6) 


where  is  the  group  delay  of  the  fixed  complex  pole.  The 
value  of  k  is  one  when  the  added  poles  are  simple.  Again  using 
Newton-Raphson  to  solve  (4.5),  the  value  of  k  can  be  Increased 
from  one  to  two  in  about  5  single  iteration  steps.  This  makes 
the  added  poles  double. 

1.  The  remaining  poles  of  the  low  pass  filter  can  now  be  Included 
in  the  group  delay  to  give 


N 


+  T 


fp 


M-2 

I. 


i=l 


2  2 
U)  +  s^ 


T 


(4.7) 


where  M  is  the  number  of  filter  poles.  The  value  of  K  can  be 
gradually  changes  from  zero  to  one  using  the  Newton-Raphson 
technique  to  solve  (4.5).  This  process  shifts  the  added  pole 
locations  such  that  the  equiripple  group  delay  is  still  main¬ 
tained. 

g.  The  conpensation  is  completed  by  noting  that  the  group  delay 
of  one  double  order  pole  at  s  =  -  is  the  same  as  one 

simple  pole  at  s  =  -  +  ju)^  together  with  one  zero  located 

at  s  =  +  J(n^.  This  latter  combination  gives  the  all-pass 

network. 

For  convenience  two  classifications  are  given  to  the  compensator 
poles.  Type  A  is  given  to  the  complex  pole  pairs  whose  imaginary  part  is 
less  than  the  cutoff  frequency.  Similarly,  Type  B  is  given  to  those  poles 
with  an  imaginary  part  greater  than  the  cutoff  frequency. 

From  the  digital  computer  program  written  to  perform  the  indicated 
steps,  two  different  all-pass  compensators  were  obtained  fromthe  Butterworth 
filters  of  degrees  4  through  8.  One  compensated  the  group  delay  up  to  the 
cutoff  frequency  (one  Type  A  section) ;  the  other  compensated  the  group  delay 
up  to  and  beyond  the  cutoff  (one  Type  A  and  one  Type  B  section) .  The  cor¬ 
responding  pole  locations  are  given  in  Table  4.  The  impulse  responses  of 
the  Butterworth  filters  (degrees  4-8)  with  (a)  no  compensation,  (b)  one  A  sec¬ 
tion  and  (c)  one  A  and  one  B  section  are  shown  in  Figures  25,  26,  and  27 
respectively.  The  Butterworth  filter  of  degree  3  with  Type  A  compensation 
will  not  work  by  this  method  because  the  dominant  pole  is  not  close  enough 
to  the  imaginary  axis  and  the  addition  of  another  complex  pole  pair  results 

in  a  delay  without  ripples.  This  particular  case  has  already  been  solved  by 
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O'Meara  with  a  single  real  pole  all-pass  section  and  needs  no  further  dis¬ 
cussion. 

Several  conclusions  which  in  turn  agree  with  the  idealized  time 

3 

domain  studies  of  Bangert  can  be  drawn  from  an  inspection  of  Figures  25-27. 
Some  of  these  conclusions  are: 

a.  Type  A  compensation  doesn.'t  necessarily  decrease  the  magnitude 
of  the  main  side  lobe,  it  merely  puts  the  large  side  lobe  in 
front  of  the  main  pulse  instead  of  behind  it. 


+1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1  +1 


79 


TABLE  4 


The  following  table  gives  the  pole-zero  locations  of  the  compensating 
all-pass  network  and  the  filter  being  compensated  (t^  =  lov/  frequency  mean  delay, 
At  =  peak  deviation  from  T^) . 


Compensator 

pole-zeros 


Filter 

poles 


Comments 


real  part 

imag.  part 

real  part 

imag.  part 

t .688654 

■^•.2.30652 

- 1 , 00000 
-  .500000 

0.0 

t. 866025 

Butterworth,  N=3,  IB  section 

T  =2.936,  At=,4611 
o  ’ 

.769347 

•t .  .420376 

-  .382683 

-  ,923790 

i. 923880 
382684 

Butterworth,  N=4,  lA  section 
Tq=6.630,  At=. 01338 

t. 587409 
t-.  386700 

■t-  .371540 
■i^l. 42347 

-  .382683 

-  .923879 

t  .923880 
t. 382684 

Butterworth,  N-4,  lA  and  IB 
section,  Tjj=8,270,  At=, 08257 

t. 500331 
t .329566 

t  .374832 
^1. 43535 

-1.00000 

-  ,309017 

-  .809017 

0.0 

i. 951056 
t. 587785 

Butterworth,  N=5,  lA  and  IB 
section,  7^=9,180,  At=.2155 

t. 644130 

t  .411718 

-1,00000 

-  .309017 

-  .809017 

0.0 

t. 951056 
i. 587785 

Butterworth,  N=5,  lA  section 
To=7.711,  At=. 06597 

t-. 546426 

t  .404039 

-  .258819 

-  .965926 

-  .707107 

t. 965926 
t. 258819 
i. 707 107 

Butterworth,  N=6,  lA  section 
Tp=8.770,  At=.1735 

t, 433629 
t. 287087 

t  .374337 
^1. 43698 

-  .258819 

-  .965926 

-  ,707107 

t. 965926 
i. 258819 
i. 707107 

Butterworth,  N=6,  lA  and  IB 
section,  T  =10.09,  At=,4094 

t .469585 

t  .397092 

- 1 . 00000 

-  .222521 

-  ,623490 

-  .900968 

0.0 

■t. 974928 
t-. 781831 
t. 433885 

Butterworth,  N=7,  lA  section 
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Butterworth  filter 


Figure  26.  Impulse  response  of  the  Butterworth  filter  with  1  A  section  of  delay  compensation 
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b.  The  time  from  the  excitation  to  the  time  the  pulse  has  decayed 
to  a  small  value  is  Independent  of  the  delay  compensation. 

c.  With  Type  A  and  B  sections,  the  main  pulse  and  the  two  large 
sidelobes  are  very  close  to  sin  x/x. 

d.  The  distortion  terms  predicted  in  the  introduction  are  of  the 
form  (assuming  lOO^^  efficiency) 

.  “  r.  m 

m  P  sin(t  -  T  +  2t  — 

distortion  =  i  (-1)  C  ,35  - ^ (4.8) 

(t  -  T  ^  2t  ^)  -' 

o  on 

Using  the  8th  degree  Butterworth  filter  with  one  A  and  one  B 
section  the  above  is  (m  =  6,  n  =  16,  T  =  11.9  and  £  =  .08) 

sin(t  -  T  +  8.9) 

distortion  =  -  .0186  - r- - ^  q  a\  •  (4.9) 

("t^T  +0,9) 

O 

Without  actually  plotting  the  sum  of  the  distortion  term  and 

.35  sin(t  -  T  )/t  -  T  ,  it  can  be  seen  that  this  will  give  a  very  close 
o  o' 

picture  of  Figure  27  (curve  n  =  8),  since  not  only  are  the  sign  and  the 
multiplying  constant  of  the  distortion  term  about  correct,  but  also  the  magni¬ 
tude  of  the  delay  shifts.  From  these  results  it  is  observed  that  the  dis¬ 
tortion  terms  so  arrange  themselves  that  the  physical  system  does  not  display 
an  anticipatory  effect. 

One  more  example  of  compensation  will  be  given.  This  is  a  5th 
order  Chebyshev  filter  (1/2  db  ripple)  compensated  with  three  A  sections  and 
one  B  section.  The  impulse  response  is  shown  in  Figure  28.  Indeed,  for  this 
example  the  output  is  very  close  to  sin  x/x  lor  a  considerable  time  on  both 
sides  of  the  main  lobe. 

Other  filters  up  to  8'th  order  have  been  compensated  by  this  method. 

One  difficulty  was  encountered  with  large  ripple  Chebyshev  filters  (2  db  ripple). 
The  group  delay  of  this  filter  is  not  a  monotonically  increasing  function  for 
frequencies  less  than  cutoff.  Several  small  local  minima  and  maxima  are  present. 
This  necessitates  an  incremental  change  in  K  of  (4.7)  to  be  very  small.  Other¬ 
wise  the  computer  program  will  evaluate  the  wrong  mlnlmax  and  hang  up.  Because 
of  this,  considerable  computer  time  is  required  (about  4  minutes  on  the  CDp  1604 
for  a  sixth  order  filter  with  5  or  6  section  of  compensation) , 


I 


Figure  28.  Impulse  response  of  a  Chebyshev  filter  n  =  5,  1/2  db  ripple,  with  3  A  and  1  B  section  of 
delay  compensation;  To  =  21.374,  At  =  .7671,  v  =  3.58< 
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Problems  of  the  preceding  type  can  be  worked  much  faster  if  inter¬ 
polation  schemes  are  used.  As  a  parameter,  K  in  (4.7)  for  example,  is  changed 
the  past  pole  locations  are  remenibered.  From  these,  linear  or  parabolic  inter¬ 
polation  can  be  used  to  estimate  future  locations.  This  method  has  been  used 
on  problems  of  the  above  type  with  very  favorable  results. 

In  summary,  this  chapter  has  made  use  of  the  digital  computer's 
ability  to  perform  the  many  iterations  required  in  going  from  a  rough  guess 
to:  a  useful  or  desirable  solution.  Indeed,  by  modifying  Example  (4.2)  slightly, 
the  equiripple  group  delay  problem  could  have  been  solved. 


86 


5.  SUMMARY,  CONCLUSIONS  AND  FURTHER  PROBLEMS 

A  rational  function  which  approximates  a  polynomial  in  an  equiripple 
manner  has  been  developed.  The  limitations  of  the  Chebyshev  rational  function 
of  Helman  and  Bennett,  mentioned  in  Chapter  1,  have  been  overcome.  Some  dif¬ 
ferences  in  these  two  functions  are:  the  old  form  has  arbitrary  pole  locations 
and  the  new  one  has  arbitrary  error  maxima  and  minima;  no  transformations  are 
involved  in  the  new  rational  function,  hence,  the  coefficients  of  the  numerator 
and  denominator  are  simple  algebraic  combinations  of  the  maxima  and  minima 
error  frequencies.  This  is  in  strong  contrast  to  Bennett's  form  which  depends 
upon  a  radical  in  a  transformation.  Consequently,  the  numerator  of  his  rational 
function  is  cluttered  with  radicals,  the  even  powers  of  which  give  a  polynomial 
and,  hopefully,  the  odd  powers  cancel. 

Several  examples  were  given  to  illustrate  the  use  and  flexibility  of 
the  new  rational  function.  Some  examples  were  used  which  have  been  treated 
using  the  old  rational  function;  namely,  the  group  delay  and  magnitude  squared 
function  approximations.  These  problems  were  adequately  handled  by  the  new 
rational  function.  In  addition,  several  problems  were  undertaken  which  were 
handled  easily  by  the  new  rational  function.  These  in  turn  would  be  difficult 
or  impossible  to  solve  with  the  old  rational  function.  Examples  of  this  latter 
type  are  given  in  the  sections  on  even  and  odd  network  function  parts  approxi¬ 
mations  . 

It  can  be  concluded  that  the  new  rational  function  can  be  a  useful 
tool  in  approximating  a  low  order  polynomial  to  give  an  equiripple  error.  This 
is  especially  true  when  less  than  the  maximum  number  of  ripple  is  desired,  or 
when  the  placement  of  the  maxima  and  minima  points  is  important. 

As  stated  in  Chapter  1,  there  are  many  problems  which  cannot  be 
adequately  handled  by  either  equiripple  error  or  rational  function.  However, 
many  of  these  problems  are  of  a  type  that  can  be  solved  using  iterative  tech¬ 
niques  with  the  aid  of  a  digital  computer.  A  numerical  scheme  has  been  suggested 
to  accomplish  this  latter  goal.  The  method  depends  upon  an  initial  guess  which 
is  already  in  itself  a  crude  solution  to  the  problem.  This  in  many  cases  is  a 
disadvantage. 

There  are  many  new  and  as  yet  unsolved  problems  which  have  been 
brought  into  focus  by  this  work.  Some  of  these  are: 
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1.  The  non-linear  simultaneous  equations  for  the  linearly  in¬ 
creasing  group  delay  need  to  be  investigated  and  solved. 

2.  In  this  work  it  was  assumed  that  various  coefficients  of  the 
numerator  could  be  set  equal  to  zero.  Necessary  and  suf¬ 
ficient  conditions  are  needed  to  say  whether  this  is  always 
possible^  or  if  not,  when. 

3.  An  extension  of  the  developed  theory  to  include  the  approxi¬ 
mation  of  piecewise  smooth  polynomials  would  be  useful. 

4.  Several  degrees  of  freedom  are  used  in  making  a  transfer 
function  realizable  by  a  symmetrical  network.  The  degree 
of  freedom  not  utilized  can  be  used  to  satisfy  other  con¬ 
straints.  One  useful  constraint  on  these  freedoms  would 
be  to  require  an  equiripple  group  delay  solution,  (however, 
the  maximum  number  of  ripples  would  not  be  present) . 
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